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FOHFViORD 


Thla la a combined report for the two projects at the Massachusetts Institute of 
Technology which are sponsored oy the Office of Naval Research under Contract N5orjb0. 


Project or. Machine Methoda of Compulation and Numerical Analysis 

This Project Is ar. outgrowth of the activities of the Institute Committee on Machine 
Methods of Computation, established In November 1950. The purpose of the Project Is (l) to 
Integrate the efforts of all the departments and groups at M.I.T. who are wording with modern 
computing machines and their applicattors, and (2; to train men In the use of these machines 
for computation and numerical analysis. 

People from several departments of the Institute are taking part in the project. 

In the Appendix will be found a list of the personnel active In thla program. 


Project Whirlwind 

Thla Project makes use or the racllltlea of the Digital Computer Laboratory. The 
principal objective of the Project la the application of an electronic digital computer of 
large capacity and very high speed (Whirlwind I) to problems In mathematics, science, engi¬ 
neering, simulation, and control. 


The Whirlwind I Computer 

Whirlwind I la of the high-speed electronic digital type, In which quantities are 
represented as discrete numbers, and complex problems are solved by the repeated uae of 
fundamental arithmetic and logical (l.e., control or selection) operations. Computations 
are executed by fractIona 1-mlcrosecond pulses In electronic circuits, of which the prlncl- 
^ ones are (l) the flip-flop, a circuit containing two vacuum tubes so connected that one 
* me or the other Is conducting, but not both: (2) the gate or coincidence circuit; (3) the 
n.agnet Ic-core memory, in which binary digits are stored as one of two directions of magnetic 
flux within ferro-magnetlc cores. 

Whirlwind I uses numbers of 16 binary digits (equivalent to about 5 decimal digits). 
This length was selected to limit the machine to a practical size, but It permits the computa¬ 
tion of many simulation problems. Calculations requiring greater number length are handled by 
the iae of multlple-length numbers. Rapid-access magnetic-core memory has a capacity of 
52,768 binary digits. Present Bpoed of the computer la iJO.OOO alngle-addresa operations per 
second, equivalent to about 20,000 multiplications per second. 


PART I 

Machine Methods of Computation and Numerical Analysis 

1 . general comments 

The 3ummer quarter marks a turning point in the 
activities of the Project. In place of the personnel who re¬ 
ported their final results at the end of the last quarter, new 
personnel are being introduced to the computing facilities of 
the Institute. As always, these new members have special in¬ 
terests in various sciences and common interest in the tech¬ 
niques of modern computing machinery. 

Moot of the reports in the following pages represent 
projects that continued in progress over the summer in the 
fields of mathematics and physics. One new and one final re¬ 
port on projects initiated in the engineering and mathematics 
departments are also included. The common objective of each 
of these is to extend the techniques of numerical analysis and 
computing into some field of application. In achieving this 
objective, the techniques themselves are enriched with useful 
sub-routines which then become available to other workers. 


!> 
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It will be recalled [l] that this problem is con¬ 
cerned with the solution of the matrix equations 


X^A (h)a = o 


V 

TlZl, the A r d6Pend UPOn the ener6y parameter h through com¬ 
binations of spheroidal wave functions. The eigenvalues h 

are determined by the usual condition that the Infinite deter¬ 
minant 




The procedure that was followed was to find the "eigenvalues" 
o a set of n x n determinants of Increasing size In the hope 
that the successive approximations h£ n ) would converge rapidlv 
as a function of n. f j 

The result was unexpected. It was found that for 
the well distortions of interest (£-.7 where £ = Interfocal 
distance/major axis) the values of h^ 

• • rCJy 


dl< 3 (i) t change by roore 


Here h ki 13 the ^'th 


than a per. cent or two for 1 # n* 4 . 
zero of the diagonal element A,, . 

KK 

Since It took better than fifteen minutes for Whirl 


wind to find a typical h£{> as contrasted with 4-5 minutes 
for the corresponding first-order approximation, I decided to 
be content with the lesser accuracy In the hundred or so 
eigenvalues which were required. 

Computations are about 80 per cent complete for a 
nucleus of A = 160, well-depth 40 mev, and nuclear radius 

x 10 cm. Five different eccentricities are being 
used, ranging from £ - .5 oblete to £ = .8 prolate. The 
energies thus obtained will be used to compute the variation 
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in the total energy as a function of the distortion. 

Jack L. Uretsky 

Reference: 

[l] Machine Methods of Computation and Numerical Analysis, 
Quarterly Progress Report No. 15, March 15 (1955) p. 19. 
Also Report No. 16, June 15 (1955) p. 21. 

ATOMIC WAVE FUNCTIONS AND ENERGIES 

The program outlined in Quarterly Progress Report 
No. 16 has been abandoned and the computation of energies has 
been completed using a lengthier but more straightforward and 
accurate procedure. 

The method mentioned in the preceding Progress Report 
gives six significant figures for the first term of non- 
lsoelectronic sequence but loses accuracy rapidly for the 
higher terms. The method used in its place consists simply 
of first holding "b" and "c" constant and finding the best 
value of "a" to within *0.02 and then doing similarly for "b" 
and "c" with A"b" = ^0.02 and A"c" = + 0 . 003 . This process 
is then iterated as many times as is necessary to give six 
significant figures. The process for one atom takes about one 
minute of machine time. The computation of the energy is 
divided up into an "a" dependent part, a "b" dependent part, 
etc. so that we need not go through the full energy evaluation 
procedure each time as we go along one of the parameter axes. 

Listed below are some results which are compared to 
previous computation done on I.B.M. W = energy in atomic units. 
(ls Z 2s Z 2p 2 ) J p 


I.B.M. Results 



a 

b 

c 

w 

C I 

3.51 

3.18 

0.965 

75.256 

N II 

3-24 

2.892 

1.027 

107.655 

°III 

3.01 

2.55 

1.003 

146.080 
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Whirlwind Results 


0 . 91*2 

1.007 

1.023 


75.257 

107.657 

146.081 


taucht The ,'’ P ° era “ ln e ° f U« variational procedure haa 
tough the writer of thla report one „aJor leaaon- „,„ely, 

t la usually beat to use the .oat stralshtforward pro¬ 
cedure possible instead or .ore elaborate ones whose value 
is evident only in special cases. 4 gl>no . , t the pr , vlous 

progress reports on this proble. will bear the above state.ent 


Arnold Tubis 


ENERGY BANDS IN GRAPHITE 


The tight-binding calculation of a two-dimensional 
graphite model is being continued. As was described in a 
previous report [i] , three-center integrals were found to be 
necessary and these are being evaluated. These integrals 
are being done by the method of expanding the orbitals around 
another center and expressing the result as a series of Gegen- 
bauer polynomials with numerical integrals as coefficients 
For several cases the series were found after ten terms to’ 
be only moderately convergent. By empirically observing that 
the coefficients in a given series are roughly in a geometric 
progression, it is possible to greatly enhance the convergence 
by approximately summing the infinite tall of the series using 
expressions derived from the generating function of the Gegen- 
bauer polynomials. To avoid numerical errors, a computer pro¬ 
gram has been written to accomplish this latter process and is 
currently being tested. 

Fernando J. Corbato^ 


9 
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Reference: 


[l] F. J. Corbato^ Quarterly Progress Report, Solid-State and 
Molecular Theory Group, M.I.T., July 15, 1955, p. 8. 

COULOMB WAVE FUNCTIONS 


The methods of getting the regular and Irregular 
Coulomb wave functions outlined by Abramowltz (Phys. Rev. £ 8 , 
1955) requires a knowledge of the Irregular function g L for 
L = 0. Up to now the problem has been one of computing the 
Irregular function 

J** 

with sufficient accuracy to be useful In any iterative or 
interpolatlve procedure for finding other Coulomb wave func¬ 
tions. To this end new subroutines were written for all the 
functions which occur In the above integrand. These routines 
are accurate to 8 and 9 places for almost all relevant values 
of the argument. One additional difficulty is the fact that 
the second term in the integrand is not well behaved as t—> 1 . 
It oscillates more and more rapidly In this limit, and the 
question arises, how to get a good estimate of its contribution 
to the whole integral. Dr. Abramowltz, when he was at M.I.T. 
during the summer for a week, suggested the following device, 
which we think largely eliminates this difficulty. From the 
second term of the integral (taken along the transition line 
f* = 2 ? for convenience, and this can be done for any value of 
^ond >} ) subtract and add the limiting form of the integrand. 

l.f. [ t — i — * E l =• 1 — i Lot} TL 

I 0 ^Uf[t - i Ley ] J- J.t = 

L irf]) - <**( z^[l «• 
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The second Integral can be done analytically; 

it 9*1 

express ion'numerlcally T ** ,bove 

oscillation., y ’ U n ° W haS the P r °P er, ty that its 

the t T ra ^ Vanl£h nCar thC 6nd P° ln t t = 1, for 

iL Z IT 3 10 thG lnt6Srand a PProach each other. Thus 
although the above integral ’ 

analytic form i* < g 1 has a somewhat more complicated 

™’ ™° h better behaved th »" the original 

integral and one should be able to 

estimate get a much better numerical 

estimate of it. T he integral .a™ can b e don. once a „d for 

The numerical technique we have In „i„a ls a aau ,. lBn 

ri:::; y 2 *? »• -* *■—-* 

.tor .^ th * -“ i "“ ~ « 


Zoltan Fried 
Aaron Temkin 
Arnold Tubls 


CRACK GROWTH IN A DUCTILE MATERIAL 


The stress pattern within a bar with an axial crack 
er torsional loading must be found as an Initial step ln 

mat! 3 , f atl n g ^ Pr ° CeSSeS by whlch tracks grow in a ductile 
r 81. Due to the complexity of the problem because of 

o/a M*h tl0n T" Unearlty ln the P laatl <= region, the use 
of a high-speed computing method is necessary. The problem is 

being studied In order that „ may be presented In a To™ 
adaptable to the Whirlwind machine. 

Joseph B. Walsh 
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THEORY OF STOCHASTIC PROCESSES 

The purpose of this report i3 to present a theorem 
which has arisen in the theory of enchained stochastic processes. 
Some ideas relating to this theorem have been presented pre¬ 
viously in these reports [4] and have found application in the 
theory of cosmic ray showers [2,3]. 

Let (A, 1 ^, P) be a probability space and let (T, Q>), 

( y be measurable spaces. The symbols D ,F" , / represent 

some general spaces made up of points ao, G, g, respectively. 

For each real x^ 0, let Qf (x), c $' (x) be measurable functions 
defined J"1 to r , y , respectively. The collections 

= {"¥(*)j o £, jl oo } will be called stochastic processes. 

In order to state the theorem of interest in this 
paper, it will be necessary to make the following assumptions: 

(1) & is a Markovian process. 

(2) The (T-field of u> events induced by *5* (x) is 
contained in that Induced by Qf (x) for all x 0. 

(3) For all x, y, x 9 0, y ?• 0, T(x) = T (y) if 
and only if & (x) = & (y). 

(4) For all x, y, x » y ? 0, & (x) = <J/(y) if and 
only if for all z, x -? z ^ y 2-0, & (x) = &(z). 

(5) The space y to which “TCx) is defined is de¬ 
numerable and g < ^for all g e y. 

(6) Corresponding to the process & there exists for 
all x, y, x > y ^ 0, a transition function, P 1 (B f , , x|G, <j), 
such that (1) as a set function of B^. £ ^^^lt is non-negative 
and completely additive, (2) as a point function of G eT it 
is measurable with respect to Q r and (3) for each set B r< f 
and each x, y, x >y >0, P^Bp, x l&(y\ y) is the conditional 
probability of the event & (x) f B p given #(<j) with probability 1. 

(1) x / Q o 0) satisfies the linear, temporally 

homogeneous diffusion equation, 

£T?(T5 r ,x/ v )= Cr 0) o) TT( G ^ r ) 

- %o)*a 
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G 0 e F 0cx < 00 , where 

^ (& ) = ^ ( ^ * 1 /( , <? / 15 r 


IT (Cr -* (r) = £> 


cxU6) = TT (& -> r ) . 


From these assumptions it follows that there is a 
transition function, P 2 (g f x^), corresponding to the process 

* ' Pl3yS the anal ytlcal counterpart of the conditional 

probability of the event T(x) = g given T(f>. The purpose 
of the theorem of this report is to establish that a certain 
condition on the function tf(Q) i s sufficient for insuring 
that P 2 (g, X |g 0 , 0) has a specified analytical form. 

It is necessary to define a special stochastic 
process. This process is simply a tool for computing the 
chance that if y arrived at the state g for some x > 0, then 
it did so by passing through a certain specified sequence of 
points of the denumerable space y . Let r(g) represent the 
apace of all finite sequences of distinct points in / ending 
with g. That is, if p* 7r(g), there is some finite sequence 
g 0' g l'**-' g N-l» g (where N may be different for different p) 
ouch that p •_ (g Q , gj,..., g). For each g ey , define 

a random variable 'P(g) such that ?(g) = p = (g , g ,..., 

g N_l» g) II and only if (l) for some set (x., x„,..., x ), 

X k > \-l > °> *?"(*) = Bq, y (x 1 ) = g 1# ... ‘Tix ) „ g ’ 

°F (* N ) = g and (2) for all x, x N > x £0, T(xj = g^ for 
some k o 0, 1,..., N, g N = g. i n words, 9(g) , p when the 
sequence of points which e T has passed through in getting to 
g is p . To complete the definition It is necessary to im¬ 
pose the normalizing condition 


£ T[ <P0i) = pJ * 1 


Now define the process 
with the statement 


1 ? 


} 
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i ( , * »«. 

is thus a random variable defined to the space of all 
pairs ( P , e), peir(a), S*/- ^TW = i» the event 

that at x V(x) = g and the process T arrived at g by passing 
through the sequence p . It is clear that for all g <X 

J p< r?c^) 

The theorem can now be stated. In the statement of 
the theorem " *" represents the operation of convolution. 
Theorem : Suppose assumption s (l)-(7) hold . If for each ^ 

cc(0) = const. = e*(e) for all &t g, then VV i£ Mar k ovian and 

0/ , , ^ f tt WCl -»f)l 


2 * P'V i- 1 J 

g N = g, where for each fixed g' (■ y ; Tf\(^'-*p >) 

is a probability measure over the space y . * 

The proof of this theorem will be contained in a paper 
yet to be published entitled, "The Concept of Enchainment--A 
Relation Between Stochastic Processes." 




Bayard Rankin 

References: 

[1] Machine Methods of Computation and Numerical Analysis, 
Quarterly Progress Report No. 14, December 15 (195*0» P* *+5* 

[2] ibid. Report No. 15, September 15 (1954), p. 48. 

[5] ibid. Report No. 14, December 15 (1954), p. 11. 

2.5 Final Reports 

CALCULATION OF NUMBERS OF STRUCTURES OF RELATIONS ON FINITE SETS 

A table of numbers of structures of dyadic relations 
has been calculated on Whirlwind I. The problem was taken up 
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primarily to test a multi-register arithmetic program for manipu¬ 
lating numbers of arbitrary length. Thus, we obtained exact 
integer answers to this problem, even though these results are 
as high as 10° . The results are given here completely written 
out, although they have primarily curiosity value. 

The problem, as described in a previous repovc^ , con¬ 
cerns dyadic relationships holding among a set of n objects. 

A complete relationship is specified by an n x n matrix of l's 
and 0's, a one in the ij place indicating that element i bears 
the relationship to element J while a zero indicates the absence 
of such a relationship. Counting the number of structures of 
relations amounts simply to counting the admissable arrays of 
l's and 0's in the Incidence matrix^ With no further restric¬ 
tions, we see that the answer is 2 n , but in this figure we 
have Included many "orbits" of isomorphic structures which can 
be permuted into one another by renumbering the objects of the 
set. The task at hand is to find how many orbits of non- 
isomorphic structures exist. Davis [l] has shown that this 
number is 

ll) str n = h Z ^ zd(7r) 

t r 

where the summand is to be evaluated for one permutation, if, from 
each conjugate class of the symmetric group of permutations on 
n objects. Every member of a conjugate class has the same 
distinct disjoint cycle scheme specified by 

(Pl» P 2 »***» P n ) 

where p^ is the number of cycles of length k in the permutation. 

The total number of conjugate classes is the number of partitions 
of n into Integral summands. The quantity b(m) is the redundancy, 
or number of member permutations in one conjugate class and is 
given by 

Pi Po Pn \-l 

b(ir) - n!(l 1 Pjl 2 P 2 l.n n P n O 
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The quantity d(ir), known as the number of "degrees of freedom" 
connected with the permutation ir, Is defined by 


n n 


d(ir) “ ZT 2Z P h p k 

h=l k=l 


graduate school research 


lrs n - XT E b ^) 2 dlrs(ir) 


d ir S ( T ) = d. 


- 2 2Z p h p k (h ' k) + k Pk 

h<.k k»l 

(h,k) = greatest common divisor of h,k 
Davis has developed other formulas for enumerating special¬ 
ized classes of relation: 

Non-Isomorphic reflexive (or lrreflexlve) relations 

ref n “ nl £ 2 dref(ir) 

? 

n 

d ref (ir) - d(v) - £ P x 

k-1 

Non-Isomorphic symmetric relations 
Sym n = nT 

ir 

d sym (lr) * E Pk {[?] + 1 + k (P k -D/2 } 
k=l 1 1 

+ E p h p k (h » k ) 

h<k 

[f]“ greatest Integer function 


Non-isomorphic antl-symmetric relations 
aSym n " HI ITb(T) 3 daaymiW) 


d asyn/ 7r ) 


it 

C Pk { + k(p k -l)/ 2 j 


+ JZ % p k < h ' k ) 

h<k 

Incidentally, note that ref n is the number of directed graphs 
n nodes and irs n is the number of non-directed graphs. 

, n J f0rmUlaS b6en evaluated ^ n ranging up 

the values are given In the accompanying tables. 

Miotic Formulae - Inspection of the various enumeration 
formulae given above shows that the dominant contribution to 
he total number of structures is due to just one of the par¬ 
titions. This partition Is the one consisting of n 1-cycles 
and corresponds to the identity transform of the group of 
transforms of the incidence matrix. Taking this term from 
each of the formulas we have 


tlons 


Nonisomorphic lrreflexlve (or reflexive) symmetric rela- 
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To show the accuracy of these approximations, we give Table VII 
as a representative table. It appears that the asymptotic 
formulae are good to about one per cent if the true structure 
number is of the order of 10 10 and are (naturally) better for 
larger structure numbers. 


M. Douglas Mcllroy 


References: 


[l ] R. L. Davis, Proc. Am. Math. Soc. 4(1953) 486 

[ 2 ] M. D. Mcllroy, Machine Methods of Computation and Numerical 
Analysis, Quarterly Progress Report No. 15 (1955) p. 10 


TABLE I Numbers of Structures of Relationships 


jail structures 

reflexive 

ref 
_ n 

symmetric 

S £V 

Trrenexlve' 

symmetric 

lrs 

n 

asymmetric 

1 

2 

1 

2 



1 — 

2 i 

10 

5 

6 

2 


2 

3“ 

104 

16 

20 

4 


7 

4 1 

3044 

218 

90 

11 


42 

5: 

2.9197 10 5 

9608 

544 

54 


582 

6 \ 

9.6929-10 7 

1.5409-10 6 

5096 

156 


21480 

7 

1.1228 10 11 

8.8203.10 8 

79264 

1044 


2.1423•10 8 

8 

4.5830 10 14 

1.7934 10 12 

2.2086 10 6 

12346 


5.7502 10 8 

9! 

6.6666 • 10^ 8 

1.3028-10 16 

1.1374 10 8 

2.7467 

10^ 

4.1594 10 11 

10 

3.4939 lO^ 5 

3.4126-10 20 

1.0926 IO 10 

1.2005 

1<T 

8.1601 10 14 

n j 

6.6603 IO 28 

3.2523 10° 

1.9564 IO 12 

1.0190 

10 9 

4.3744 10 18 

12 

4.6557•IO 34 1.1367 IO 51 

6.5234 • IO 14 

1.6509 

10 11 

6.4540 10 22 

13 

9.0169 io 40 

1.4669-lO^ 7 

4.0540 10 iY 

5.0502 

10 13 

2.6378 IO 27 

14 

1.1521 10 4 °:7.0316 10 43 

4.7057 io 20 ,2.9054 

10 16 

3.OO37 10 52 : 

15 

4.1233 IO 55 

1.2583 10 51 

1.0231 IO 24 

VO 
CM 
-=*- 
r —1 

K\ 

10 19 

9.5773 io 37 

16 

5.5343 io 6? 

8.4446 -IO 58 

4.1788 IO 27 

6.4001 

10 22 

8.5888 IO 43 
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TABLE II Numbers of Structures of Dyadic Relations 


n 

1 

2 

3 

4 

5 

6 

7 

8 

9 6666 

10 

90545 

11 

85078 

12 

89066 

13 901685 

19117 

14 

57604 

45256 

15 4 
79295 
20975 

16 

25727 

57553 

89642 


11 

45829 

62157 

49349 

18075 

46557 

03112 

91267 

62528 

74157 

77230 

•12334 

18834 

59863 

62971 

09620 

93721 


2 

969 

22829 

71000 

21539 

98391 

6660 

85386 

45648 

66511 

11300 

96061 

1152 

55389 

31424 

41401 

69376 

65854 

20722 

02145 

27245 


str 

n 

2 

10 

104 

3044 

91968 

28992 

08928 

61728 

27936 

3493 

61856 

34219 

36288 

25869 

04256 

76041 

48096 

05015 

34617 

28672 

68606 

48648 

35136 

5534 

05192 

19348 

80352 
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TABLE III Numbers of Structures of Reflexive (or 
irreflexive) Dyadic Relations 


n 




ref n 

1 




1 

2 




3 

3 




16 

4 




218 

5 




9608 

6 



15 

40944 

7 



8820 

33440 

8 


179 

33591 

92848 

9 

13 

02795 

68243 

99552 

10 

341260 

43195 

29725 

80352 

11 



3 

25229 


09385 

05588 

61111 

97440 

12 


11 

36674 

54308 


25400 

57443 

38940 

04224 

13 

146 

69085 

69271 

29298 


69037 

09607 

53162 

20928 

14 

56566 

15234 

99952 

7031 

13855 


06555 

97990 

40912 

17920 

15 


12 

58345 

26155 


04488 

67281 

04228 

58105 


99188 

12349 

03206 

83OO8 

16 

8444 

60738 

34225 

80541 


87807 

17815 

32315 

89171 


86915 

O3432 

37883 

67872 
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table IV 


Numbers of Structure 
Relations 


s of Symmetric Dyadic 


n 

1 

2 

3 

4 

5 

6 

7 

8 

9 


10 



11 



12 



13 


405 

14 

4 

70568 

15 

IO230 

63423 

16 

41788492 

03082 


sym n 

2 

6 

20 

90 

544 

5096 

79264 



22 

08612 


1137 

43760 

1 

09262 

27136 

195 

63634 

35360 

65233 

50845 

92096 

40227 

34209 

96800 

64216 

11199 

63904 

47118 

94310 

54720 

O2323 

60582 

29792 
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TABLE V Numbers of Structures of Irreflexive (or 
reflexive) Symmetric Dyadic Relations 


n 

1 

2 

3 

4 

5 

6 

7 

8 
9 


10 


11 


12 


13 


14 

29 

15 

31426 

16 

640 01015 


irs n 

1 

2 

4 



2 

11 

34 

156 

1044 

12346 

74668 


120 

05168 


10189 

97864 

16 

50911 

72592 

5050 

20313 

67952 

05415 

56572 

35488 

48596 

98043 

08768 

70452 

75578 

94928 
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TABLE VI Dyadic S Relations^ UreS ° f A " tls »ric 


9 


41 

10 


81600 

11 

4374 

40620 

12 

39836 

93872 

13 

35571 

22500 

14 


300 


80530 

05349 

15 

957 

72686 


49990 

83757 

16 

84182 

49161 


38402 

27902 


21 

5750 

59392 

74490 

99707 

07497 

263 

90533 

36589 

84908 

34898 

92075 

16546 

32471 


asym 

n 

1 

2 

7 

42 

582 

21480 

42288 

16219 

43052 

11040 

47314 

645 

39356 

77967 

73136 

61589 

93399 

11549 

81003 

8588 

12893 

44414 
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TABLE VII Comparison of Asymptotic Structure 
Formulae with True Formulae 


approx. 
_ value 

str 2 1.117 10 1 


7 

. true 
va lue 

1.123 io ] 


approx. 

value 

3.493 10 * 
3.411 10 Z 
■993 10 1 
.970 10 7 
8.140 10 1 


10 

n = 

15 

true 

approx. 

true 

value 

value 

value 

3.494 10 2? 

4.123 10 55 

4.123 10 

3.413 10 20 

1.258 10 51 

1.258 10 

1.093 10 10 

1.016 10 24 

1.023 10 

1.201 10 7 

3.102 10 19 

3.143 10 

8.160 10 l4 

9-577 10 37 

9.577 io : 


Project Whirlwind 


1- REVIEW AND PROBLEM INDEX 


During ?l7 1 , d h 0r Jun * 13 ' t0 September 18, 1955 

computer time allocated to the Sclentif?r h ? urs ?' the 562.2 hours of Whlrlilnd I 

remaining 22U. 9 hours of the allocated Com Putationa (S&EC) Group The 

call-ration, demonstrations. ?ape conve^s^ns ?or U f? d ‘? r ? emlnsl equipment «£“ 

run operations not logged to apecinc pro^ems " Uboralor ^ «"d various Inte?- 

problem. !» f‘ elds 0f a PPl‘c«U°n». The results of 18 of 

toral theses. 1 .^Wrlng^ ^ster's^and^cn 1 ^ 0 ”' 0f ^L^'r^entdo^ 
8.ve originated from research 

nave been submitted, ^n th<Tn ra^tah^-Lhe* to ,1 he Problems for which progress reports 
th appI J? atlon - and the source and the amount P n? are arranged according to the field 

academic credit and S ^ 
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PKOBLEM INDEX 


iuil«*l 

lug 

Aet-oeiaallc and 
Structures Research 
Utontorjr 


Inn. Aoalni stratum 
Chemical Engineering 


Civil Engineering 


Oeolugy *nd 
Oeopnyatca 


Lincoln uMratun 


4® rvomecitan l an 
Laboratory 

Miscellaneous 


pa script urn 

Transient response of aircraft atruiiure* to aerodynamic haatln« 
Hurlumta* aUMUaar woes. ahapa*. frequence* 

Tranalvnt temperature of a ooa type «■#*• 

Factoring high order polynomial* 

Traiextory calculation* for • recast during powered riltfil 

* fart tally continuous wooden teams 

Nuswrically vontrollad all Uni hcMw turbine bi#d* 

Tranalanta In distillation columns 
Calculations for tnv HIT reactor 
Haat tranafar In turbulent rioa 

Critics. aiaaa calculations for cylindrical geometry 
Prediction of cnrometograpnia aap*r*tlona 

t»na*l« buckling 

Dyrvaa.c analysis of an aircraft Interceptor 


Solitary wav* generaling can 


Outdance and control 

Oat* reduction for 1-1 fir* control 

.lyatem analyata 

Klganvalu* problem for propagation of electromagnetic wave* 

ionosphere confutation 

Ueneral nayd tat solution 

Tropoapnerlc propagation 

Error analysis 

Nadar correlation 

Plow of coagiroaaltil* fluid* laarothenaopreasor) 

Laminar boundary layer of a steady, compressible flow In tr.* 
entrance region of a tube 

3ynoptto cllaiatalogy 

•Towputatlona of tna fleida or vertloai velocity and horlaontal 
divergence 

investigation of in® vortleity rield In tr.a genera, circulation 
of the atanspner* 

Correlation function 

Spertr*. analysis of aiasiperlt data 


Self-consistent aolecular orbital 

An lugentet plane save ea’.nol as spyltad to sodlus 
Exchange Integrals uataeen real Slater orbital* 

Tran afore**, ton of tntagra'.s for diatomic soleculea 

Neutron-deuteron scattering 

Atomic Integrals 

Theory of neutron reaettuna 

Application Of the AM* method to Body- and faca-cenlerad iro 
Energy tsveia of diatomic hydrides 
Evaluation of teo-center aolecular tntegrala 
Electron dlffuelon in an a 1 #ctroawgne 11 c field 
Analysis of air anoeer data 

•kilttpl* acattarlng of eaves from a spatial array of 
apnerlcat scatter*™ 

Energy levels or d.atomic hydride* UH 

Augmentad plane eave method aa applied to crystal cnroeilua 
Atomic eave function* 

Electronic energy of the nellum aoievular Ion 

Nelatlvlatle atomic eavefuncHone 

Ture and impure AC* cry alal 

tmta reduction program, polynomial fitting 

Subrout lnea for the nuMrlcally control led milling maonir.e 

comprefwnalve ayetem of service routine* 

S h EC Subroutine study 
Ultrasonic delay lines 

irlaon of simplex and relaxation methoda in linear 


rrctlee 

Min of 

or 

Number 

7*1 Time 

frograaver 

e36 C. 

368.7 

L. Schmidt 

7H C. 

3V.7 

E We law re 

179 C. 

*9.6 

L. Schmidt 

11* C. 

20.6 

V. Howard 

310 e. 

60 .> 

J. /rlgg* 

302 

11.B 

N Weinstein 

267 a. 

190.5 

0 Bromflald 

2*1 B.R 

i®:? 

a. Davla 

266 A. 

R Trooat 

29-1 C, 

373.4 

A. Turano 

3 ‘0 8. 

•99.9 

j Tow* 1 1 

30) C. 

•7 **.§ 

J. Fischer 

791 N 

7*.e 

R. Jones 

75<J C. 

174,5 

K. Mvanagh 

264 C. 

7iS :l 

J. farnnla 

7/1 B 

F Engel 

106 C. 

767.0 

S. Simpson 

761 C. 

377. • 

N. (kiergwr 

Ill N. 

127.1 

J. HtMAsley 

2J9 C. 

234.9 

J. Unlng 

2*4 C. 

7.4 

J. Stare 

2‘7> C. 

ic;.i 

R. M>nl 

IV) l. 

723.2 

H Dwight 

259 L. 

304'..? 

t> Brennan 

373 L. 

426. • 

g. Rotanberg 

30b L. 

282.1 

H. Uelgnt 

312 L. 

255.8 

l. Shapiro 

316 L. 

17.0 

W kelnat*In 

121) B.R. 

or .6 

A. grtekaen 

199 ■ 

110.3 

T. Toong 

155 N. 

2*6.9 

E. Kelley 

72* N. 

40). 3 

J. Austin 

726 D. 

382.6 

R. Ffafrer 

28 u t 

•tf.6 

t. Hanna 

306 D. 

43.9 

8. Saltzman 

29/ 8 . 

371.4 

J. Baron 

30/ C. 

39.1 

J. Baron 

1*4 N. 

60.9 

R. Neabet 

19* B.R. 

277.6 

N. Saffren 

20* S. 

14.7 

r. Rerryman or 


61.1 

Untv, of Chicago 

21o R. 

8. Neabet 


225 B.k. 
23 - N 
2*5 R 
253 N. 
2tO N 
2tw R. 
265 L. 

2? 3 h 

27 a K 
770 M. 

?<* A. 

S* »’,N. 
.76 C. 

132 t>. 

100 . 

INI. 

216 C. 


Whirlwind l - Ah A ilO) tranalatlon program 
Course 6.r 

frequency analyale or aperiodic functions 
Noutlnas for course 6.601 


m.k 
22.4 
1062.8 
161.0 

218.7 

376.2 
01.6 

113.2 

81 .0 

411.9 

606.7 

246.6 

1054.4 

3$.5 

918.3 
74.3 

2584.2 

T!:I 

266.1 

*5:1 


8. Aghajanlan 

D. Arden 

0. Clark 

M. Karakaahlan 
0. Roster 
M. Saffren 
H. Neabet 

B. Molselwltach 

C. Schwartz 
L. Howland 
0. Rosa 
J.M.Runyon 

f. Helwlg 
3 A EC Oroup 
R. Bishop 

E. Eairra 

3. frankovlah 
U. Kcclea 
J. Roaeamn 
A. Siegel 


Table 2-1 Current Problems Arranged According to Field of Application 
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PHOBLKX 


Msthemmtlcal Frobism 

1. Matrix a.get ra and equations 

t-iJ a* 1100 * 01 Mona 11 gation 

Orthonormau*at Ion 
Hoot of a determinants! aquation 
Linear aquations 
Simultaneous equations 
Eigenvslues 
Eigenvalues 
Eigenvalues 
Invars ion 

2. Ordinary different;aaquations 

Seven nonlinear first order 

System 

or ' J «’ r equations 
3«t of nonlinear first order 
wave equation 
Set of eight Tlrat order 
Second order 

System or ordinary differentia; equations 
Nonlinear differential equations 
Nonlinear dlfrarentiai equations 
Second order 

Second order nonlinear dlffarentlal equations 

3 Partial differential equations 
Second order parabolic 
First order system 
Flrat order system 
Second order 

Simultaneous partial differential equations 
Second order nyperooltc partial 

4. Integration 

Integra; evaluation 

Stationary point of a variational 

Overlap integrals 

Integration 

Integral transformation 

Overlap Integral 

Integration 

Integration 

Integrals 

Freanal Integral 


**. Statistics 

Multiple lime series 

-slculatlon of the coefficients of a a 

regression system 
Multiple time series 

Maximum likelihood estimation 
Multiple tloK aeries 

6. Transcendental equations • 

Nonlinear equations 
Nonlinear equations 
Curve fitting 

- Data reduction 
Data reduction 
Surface rittlng 
Mias distribution 

8. Oroup theory 

denerat Ion of projection operators 

9. Complex algebra 

Complex roots and function evaluation 

1C. Fourier series 

Fourier synthesis 
Summing series 
Fourlar synthesis 

11. Linear prograsmlng 

Linear programing 

Linear prograsmlng 
Linear programing 

12. Factorlastion or Polynomials 

13. Algebraic equations 


ultiple 


14'. N. 
218 N 
266 A 
270 ». 

m c. 

?7B N 
20‘, N. 
2 nd N. 
312 L. 


haaic operations 
Schmidt process 
Iteration 
Crout‘a method 
Herat lor 
liiagnna I i rat l*r 
DlagonaUtailon 
diagonallsstlun 
Crout s method 


Fourth order Nutta-0111 
0111 a method 
Ql 11 a method 

".econo order Kutta-Olll <?H 1 

MlJne predict ion-corrector r,.rmu* 2*»‘ 
Olira method ,»yj 

Oauaa-Jackson method * 0 *, 

Difference equations 291 

Dili's ms trod 

Olil'■ method jy 7 

0111'a method j;o 

0111a method ?65 


170 

a? 


b.n 


Explicit finite diffarences 
Finite differences 
Finite dirferencss 
Finite differences 
Finite dirferences 
Finite differences 


Oauasian quadratic 
Simpson a rule 

Evaluation of analytic Toms 
Simpson's rule 
Algebraic recursion lormula 
Evaluation of analytic forma 
Trapezoidal rule 
Oauaa quadrature 
Barnett a Cou,ion expansion m 
Conversion power series with 
complex agreement 


236 C. 
224 N . 
226 D. 
299 C. 
393 B. 
179 C. 


20 - 1 
775 I 
760 J 
2/8 I 
7J4 > 
262 | 
304 i 

312 1 
at nod 309 j 


Prediction t>y linear opmratore H*> c. 

Inner products g. 

Correlations, means, frequency 

distribution 299 L. 

Means, weighted sums, etc. 2?1 B. 

Correlation and tranafonas 2b0 U. 


Steapest descent 264 c. 
Iteration 272 L. 
Least squares 273 N. 


Polynomial rittlng. etc. 126 C, 
Polynomial fitting, etc. 267 B. 
Arithmetic operations 2<*4 c. 


Machine generation 
Herat inn 


Direct evaluation 261 C. 
Direct evaluation 274 N. 
Direct evaluation 306 D. 


dimples method and relaxation 

method 219 

Simplex method and relaxation 

method 316 t- 

simplex method ana relaxation 

method 216 c. 

Mltphcocr's method jjq C. 

Newton-Naphaon j 07 c. 


Table 2-II Current Problems Arranged According to the Mathematlca Involved 


> 
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2. WHIRLWIND CODING AND APPLICATIONS 


2.1 Introduction 

Progress reports as submitted by the various programmers are presented In numeri¬ 
cal order In Section 2.2. Since this summary report presents the combined efforts of DIC 
Projects and 6913, reports on problems undertaken by members of the Machine Methods 

or Computation (MMC) Group have been omitted from Section 2.2 of Part II to avoid duplica¬ 
tion of Part 1. For reference purposes, a list of the MMC Group problems appears on 
page 68. 

Letters have been added to the problem numbers to Indicate whether the problem 
is for academic credit ana whether It Is sponsored. The letters have the following sig¬ 
nificance: 

A Implies the problem is NOT for academic credit. Is UNsponsored. 

B Implies the problem IS for academic credit. Is UNsponsored. 

C implies the problem is NOT for academic credit, IS sponsored. 

D Implies the problem IS for academic credit, sponsored. 

N Implies the problem Is sponsored by the Office of Naval Research. 

L Implies the problem Is sponsored by Lincoln Laboratory. 

The absence of a letter indicates that the problem originated within the S and EC Group. 




?.? Problems Being Solved 
j,; - COMPREHENSIVE SYSTEM OF SERVICE ROUTINES 


whirlwind codinq and applications 


.ormg the of e the e woSklSg S ^m«. R * POrt a “ — 

structure or^he'comprehenslvrsiislem?* 5 0 " * maniJal " hleh wlU deacrlte i n some detail the 

F. C. Helwlg 

Digital Computer Laboratory 

.u6 C. MIT SEISMIC PROJECT 

Wlon of the use of statlatH-al L anai'ysrs P technirue tl t m 106 U concerni!d with the ln- 
;.nd In particular with the separation of >efle^!ons 9 frnm^"'! 8 " 10 recura Interpretations, 
th-s.- recorcs. More complete descriptions of the ™ / background Interference on 
;; :ne V: i he Dleltal computer Jbor«o?y B?week!v P Rennr? ‘TV*’ 6 a ?P«»«hes uaed are con- 
tien of Reflections on Seismic Records by Linear J ?[} e 15 ■ ly '3. and In Detec- 

and Hurley--GEOPHYSICS, Vol. 18. No. 3, July 1 «3)f * ‘ (k,adai,orth • Hoblnson, Bryan, 

No. LI) for our recent approaches to the problem^' Durln^th 1 *]' 1 * report (Summary Report 
centrated on the following problem- the 8 ll ? e last l ’ uarter we have con- 

mo.-ram Tor possible reflection signals where ?,„!?}}* 8P r ead muUUra « sels- 

■race may vary due to both dip or reflector and ^"ttlon of the signal on each 

... o, 

increment of tl^yTrosi-co?relauCns e «ro P a d the r recoM PU Il° ae avera 8 ea - for each small 
eapect a reflection to Tall for various trial ! along eurvea °n which we would 

function. (Actually It averagesallong tio stra! nr' accor ^^ l o any given velocity 
traces, one on the lower hair, which represent tf£\.hnUL^segments, one on the upper hair 
ire.) For each Increment the program then fin's that ai S ° f actual reflection curve 
tverage. prints out this angle and the - th * ? lp sngle which maximizes this 

to an elevation and prints this elevation. ’ onverts the cenler time of the Increment 

different v^ue^orthe'aver^ng-fen^h^nd^Mi? T. ^ ™ actual ™ord. with 
suits are very encouraging. g "' and trial step-out spread parameters. The re- 

reproducUons ! ‘my C oe r obtIlned r throCgh t the special'collectlo 0 ^? 1 '^?''*' 1 " 8 compar, l'a. but 
Library, six months after the repo^ Irl TnVH ““compiles" 1810 '’ th<? Ha * da " 

~ *™ ■■>'2550 JXfW. rTLSTJJ »■ —. 

S. M. Slmpaon, Jr. 

Oeology and Qeophyalca 


l2 ° 8>N ' oJs EFPECTS 0F WAT£R injection into hioh-temperature. 

component. ^Hefan^aerothe^rietri^ d ? valo P ment a Potential g.a turbine 
hot gas stream Is brought abo™by cooling Sf^he stream *? a ^ gnatlon Pleasure of a 
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WHIRLWIND CODING AND APPLICATIONS 

discharging It to the atmosphere. At the entrance of the duct, special Injectors deliver 
minute Jets of water which are in turn atomized by the rapidly moving gas stream. 

TTie changes In state within the aerothermopressor are brought about by the si¬ 
multaneous thermodynamic and dynamic effects of (a) evaporation of the llauld water, (b) 
momentum and energy interactions between the phases, (cj friction, and (d) variations In 
the cros 3 -sectional area of the duct. Under proper circumstances, these effects bring 
about a net rise In stagnation pressure across the device. Further descriptions of this 
device may be found in earlier reports, beginning with Summary Report No. 32, Pourth 
Quarter, 1952. 

TTie role of Whirlwind I in the successful development of the aerothermopressor 
Is intimately connected with the determination of performance characteristics of the de¬ 
vice, under all conditions of operation, by means of a comprehensive one-dimensional analy¬ 
sis of the process. This analysis Involves the simultaneous solution of seven, non-linear, 
first-order differential equations. 

During the past quarter a number of computations were made using the new Whirl¬ 
wind program for the aerothermopressor (described In the Quarterly Report of March 2k, 

1955. No. 111). The computations can be divided into three groups: 

1. Theoretical calculations to compare with experimental results from the 
aerothermopressor being tested In the Gas Turbine Laboratory. These comparisons are In¬ 
cluded In an Sc. D. thesis by Arthur A. Fowle entitled An Experimental Investigation of 
an Aerothermopressor Having a Gas Flow Capacity of 25 Pounds Per Second". A typical 
example of such a comparison (Figure 19 from the thesis) Is Included with this report. 

2. A theoretical investigation of the desirability of having supersonic gas 
velocities at the plane of water Injection. Calculations Indicate performance of the 
aerothermopressor improves when supersonic velocities are used. This result, however, does 
not include the effect of increased friction upstream of the water Injection plane. 

3. A study of conical diffusers with different angles of divergence In order 
to de-iign a better experimental aerothermopressor for testing in the Gas Turbine 
Laboratory. 

The aerothermopressor program Is being carried out at M.I.T. under the sponsor¬ 
ship of the Office of Naval Research and is being directed by Professor Ascher H. Shapiro 
or the Department of Mechanical Engineering. The theoretical aspects or the problem treat¬ 
ed by Whirlwind I were programmed by Bruce D. Gavril and are now being carried out by Alve 
J. Erickson. 

A. J. Erickson 
Mechanical Engineering 
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WHIRLWIND CODING AND APPLICATIONS 


126 D. DATA REDUCTION 

Problem 126 la a very large data-reduction program for use In the Servomechan¬ 
isms LaDoratory. The overall problem is composed of many component sections which have 
been developed separately and are now being combined into complete prototype programs. 
Descriptions of the various component sections have appeared In past quarterly reports. 
After the development and testing of the prototype Whirlwind programs is completed, the 
rant Mill be re-coded for other, commercially available, large scale computers, 
(procaoly the ERA 1103» IBM 701 and IBM 704 computers), for U3e by Interested agencies 
for actual data reduction at other locations. The programs are currently being developed 
by Duuglas T. Rosa and David P.McAvinn with the assistance of Miss Dorothy A. Hamilton, 
Servomechanisms Laboratory staff members. This work is sponsored by the Air Force Arma¬ 
ment Laboratory through DIC Project 7136. 


The nature of the problem requires extreme automatlclty and efficiency in the 
actual running of the program, but also requires the presence of human operators in the 
computation loop for the purpose of decision making and program modifications. For this 
reason extensive use is made of output oscilloscopes 30 that the computer can communicate 
wltn the human, and manual intervention registers so that the human can communicate with 
the computer in terms of broad ideas, while the computer is running, and have the computer 
program translate these Ideas into the detailed steps necessary for program modification 
to conform to the human operator's decision. The program which does this translation and 
mod 1 flea’lon is called the Manual Intervention Program (MIV). The most recent version of 
the prototype data-reduction program In colled the Basic Evaluation Program. 

Moot of the computer time during this quarter was used In testing and modifying 
the WWI-1103 Translation Program described in Problem 266. Three programs were written In 
vhe ERA 1103 code for translation by this program. Two of these programs were written to 
teat the oscilloscope display unit attached to the Air Force Armament Center 1103 Computer 
at Egllri Air Force Base, Florida. One is a .'Jcope Plexowrlter Program which writes charac¬ 
ters on th*.- scope corresponding to a Plexo tape read into the computer. The other program 
is the Mouse Program which simulates a mouse solving a 16 x 16 maze with the successive 
stepa drawn on the scope by the program. The third program, A Basic Input Translation 
Irogram lor the ERA 1103 Computer, is a companion to the program of Problem 266 except 
that symbolic addresses are not allowed. This program will be used to make corrections 
to programs translated on WWI at the 1103 site. 


m tvaaic evaluation rrogram has again been rewritten to Incorporate still 
newer versions of t:he equations and to simplify the organization of storage locations. 

IhSi r?! 1 J ? rog r am haa also been a *> aln rewritten using many preset parameters to allow 
additional Inputs and outputs In arbitrary formats. The sine and cosine routines have 
Si^~- e L" 8dlfle ? 1° U8e local TayJor expansions to greatly Increase the speed of the 
SJirrTJT.Si !!? '“J* ln accuracy. The program chooses the Taylor Expansion IT It will give 
sufficient accuracy, but chooses the library subroutine otherwise, and U3es the results to 
generate a new Taylor expansion about the new point for use on succeeding steps. 

. The Manua: Intervention Program (MIV) has gone through several new revisions In 
efforts t° generate a system which will alliw, eaaentlally, proira^ng wl?h progr^ma 

Problem 100°ao U that 1 th lr ’ oporal t ton - The Dlr "‘or Tape Program has been modified under 
understooc^form?^ 8CU ° n Pa " b * ^nU^af J^JSSST' 


D. T. Ross 

Servomechanisms Laboratory 
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D - 3UBR0OTINES P °R THE NUMERICALLY CONTROLLED MILLINO MACHINE 

"" following work was done during the past garter under Problem Ija, 

rh^ri c ^;?? n? ^ e ^“^ c iSchfL per an 

::: .... 

........ snE- P syKr““ 

Problem 1 3 P la now terminated with regard to Whirlwind I. 

T. Nagle 

Servomechanisms Laboratory 

141 SUBROUTINE STUDY OP THE SCIENTIFIC AND ENGINEERS COMPUTATIONS (3 A «, OHOUP 
J ’ Factorization of Polynomials 

numerical factorization methods. 1 d 1 DCL - J '*-1- Research Is continuing on 


3 ' Iterative Solutions 

systems. The P methSd"??plrti^Ila^y sulted^o 10 ^ f ° r lteratlve solution of linear 

”‘ n ssj - 

the^matrl 8 *' ^°^«™regL?errne!', ; l ne a p roild l e.j n to a -.xrca^^^fh 8 ' In the rirat place - 

frois^i 1 w ereaa reduction processes such as Oauss JnrrtL f *5 e nor '- :! cro elementa ln 
from the character of the matrix. Second v lnt? ^ °‘ Crout cari not a ® prorit 

order ^t^overcome ^round-of f^^ne^an expect ^th* 4 ”!! omber^o^superf'tuous 8 ^ 8 ^!tfl^ln 4 * 1 

convergence for large system, 23 SK/SKi SStfi * ££j?" * 

ln r cor n raa " r lcrs e ^" Pu t C ’ n ear(he h llmning e sl« 8 (28 t x C 23ro? ravoPabJ >’ with Grout's method 

° method takes a time proportional to n3* k ?n nrol?* r0ughly Proportional to n^ while 
symmetrt 18 \ ln } ted onl r b > tbe oumber jf o?f?diagohaf ereSenf: <, ? lr fr nla - tht Iterative' 

all the characterlstlcs l or a ™°?ices of°the “ thod "Rich beat takes advantage of 

routine arc described In »o™ , d r «."l°fn t ^^d:m 1 ^L K !Sl? 0,,I ' r0nted ' "* ™th^fLthe 

M. D. Mcllroy 
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з. Histogram Plotter 

registers. M D Mc i lroy 

и. KU‘.A RAPID SIWE-COSINE 


rrs Srss ™ .xsss 

temporary storage. 


M. D. Mcilroy 


034 SCOFE MRA DECIMAL FORMAT 


This routine displays normalized generalized decimal numbers 1 rom the Jh'-e' - - 
pretei Multiple Register Accumulator (MRA) on the scope In any rectangular array that does 
not overreach a basic Uo * AO grid on the scope Race. Positioning of numbers proceeds 
seoui Mlaily eltner uy row or pv column. The Tormat pattern may be changed at will by the 
programmer by*meana of program parameters. This routine occupies 248 registers and uses 
2 registers of temporary storage. 


6. OS', SCOPE MRA DECIMAL OUTPUT 


M. D. Mcilroy 


Basically, this routine Is Identical to 0S4, except that the format Is fixed by 
preset parameters. 0S5 occupies 205 registers and uses 3 registers of temporary storage. 

M. D. Mcilroy 


7. Printout Subroutines 


The following subroutines have been added to the subroutine library: 


a. 

OD 

6 

b. 

on 

7 

c. 

OD 

8 

d. 

OD 

9 

e. 

OD 

10 

f. 

OD 

11 

«• 

OD 

12 

h. 

OD 

13 


Delayed Single Length Decimal Integer Print 
Delayed Octal Print 

Single Register Block Print (program parameters) 
Single Register Block print (preset parameters) 
(30-J.J) Decimal CS Print 
(30,15) Decimal Print (CS or ww) 

CS (30-J.J) Block Print (preset parameters) 

CS (30-J.J) Block Print (program parameters) 


Pour of the delayed output routines (OD 6, OD 7. OD 10, and OD 11) print out 
numbers. The first three are rewrites of OD 3, OD 5. and OD 2, respectively. In each 
case they are shorter than the original. All four routines have been written so that 
they may be used with the remaining four (OD 8, OD 9, OD 12, and OD 13) for printing out 
blocks of numbers from high-speed storage without magnetic tape stoppage between numbers. 
The number printout routines have banks of entrance blocks providing various terminating 
characters. These banks lie at the end of the routines. Unused entrance blocks can be 
deleted, thus shortening the routine. 


Two of the four block print routines provide for printing single register num¬ 
bers. The other two are for double register numbers. The parameters for the block 
routines are: 


(l) The Initial hlgh-apeed address of the block of numbers to be printed 
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(2) The addreas of the subroutine printing the number (0D 6, 7, 10, cr 11). 

(3) The number of numbers In the block. 

(*0 The number of numbers per line. 

(5) T*»e terminating character between numbers on a line. 

Of each pair of block routines, one has this information In preset parameter 
form, while the other has It In program parameter form. The preset parameter routines are 
snorter but less flexible. 

0D 11 Is modeled after the PA PM printout routine. It provides more correct 
-lULts than any previous LSR output. For numbers with exponents of large magnitude It la 
the most rapid routine to date and It will successfully print numbers with t Inary ex¬ 
ponents up to ±2*^. 


R. Watson 


FU RAPID SQUARE ROOT 


PU 8 is a square root routine which takes, as Its first approximation, the first 
three terms of the Tschebyacheff series and then does one Iteration by Newton's method. 

This routine takes 5 milliseconds as against 30 Tor FU 2 . and It is one octal digit less 

accurate. 


R. Watson 


1*4** N. SELF-CONSISTENT MOLECULAR ORBITAL 

Several revisions were made In subroutines developed under this problem In 
connection with the work on Problem 28b. 

R. K. Nesbet 

Solid State and Molecular Theory 
Group 


lb' N. SYNOPTIC CLIMATOLOGY 

During the past quarterly period, two phases of the Synoptic Climatology 
Project have utliLzed Whirlwind I. The method, essentially one of multiple linear 
regression, has been developed and Improved over the past two years, mainly by Dr. Thomas 
F. Malone and Mr. Robert 0. Miller. 

The first phase has dealt with the prediction of 2*4 hour prec IpJ tat Ion amount 
and duration at Boston. Massachusetts. The analysis ha3 been extended to Include the 
effects of the upper level circulation and moisture distribution, In addition to the aea- 
level pressure pattern on the variability of rainfall. 

In this study two basic operations are involved: 1) representation of the 
fields of certain continuous variables (e.g. measure, moi3lure) in terms of series of 
coefficients of orthogonal polynomials, and 2; matrix solution and inversion. Programs 
have been written previously by Mr. Robert 0. Miller to accomplish both of these opera¬ 
tions and only slight modifications were necessary to fit them to the present study. 

As an extension of the above, It la planned to apply essentially the same 
methods to the specification and prediction of 5-day mean precipitation over the entire 
United States. 

The second phase has dealt with the specification, prediction, and verification 
of forecasts of 5-day mean temperature anomalies at Atlanta, Ga.; Blsmark, N. D.; 
Indianapolis, Ind.; Little Rock, Ark.; and Portland, Me. Also, the apeclficatlon and 
prediction of 5-day precipitation totals at the same stations has been handled In this 
study. The prediction parameters used are standardized correlation coefficients between 
700 mb. 5-day mean height anomalies and the elements or Tschebyacheff orthogonal poly¬ 
nomials. This phase of the project la nearing completion and it is hoped to combine all 
of these results In a scientific report. 
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irfHINbWIND CODING AND ABDICATIONS 


Elizabeth 

Houghton, 


,^WS5‘A.*i‘« isrJM.W. 

and Dr. Edward N. borenz, Department 01 Meteorology. 


W. D. Sellers 
Meteorology 


1, . C. TRANSIENT TEMPERATURE OP A BOX-TYPE BEAM 

T „ Henort No 3* a urogram for computing the transient temperature and 

mrnimmm&Mm 

ne correlated with experimental worn recently conducted by Mr. S. J. Engel of th. 
elastic and Structures Research Laboratory. 


Luc ten A. Schmlt 
Aeroelastlc and Structures 
Research Laboratory 


lv3 L. EIGENVALUE PROBLEM FOR PROPAGATION OP ELECTROMAGNETIC WAVES 

The previous report on this problem applies at the present time. A set of 
practical results which can be compared with the experimental results will be obtained in 
the near future. 


H. B. Dwight 
Lincoln Laboratory 


B.N. AN AUGMENTED PLANE WAVE METHOD AS APPLIED TO SODIUM 

Mr. M. Saffran of the Solid State and Molecular Theory Group continues to test 
the programs mentioned in the previous Summary Report. 

It has been found that the matrix element generation routine which has been 
used Is overly slow and largely inaccurate; therefore, a new routine is being written. 

M. M. Saffren 

Solid State and Molecular Theory 
Group 


199 N. LAMINAR BOUNDARY LAYER OF A STEADY, COMPRESSIBLE FLOW IN THE ENTRANCE REGION OF A 
TUBE 

In connection with the research on heat transfer to a stream of air flowing at 
supersonic speeds in a round tube, a theoretical investigation of the characteristics of 
the laminar boundary layer in the entrance region of the tube has been carried out. The 
boundary layer equations of continuity, momentum and energy are to be solved for specific 
entrance Mach numbers and thermal conditions at the tube wall. 

Gill's method is used in the numerical solution of these equations. 

Solutions of the first three seta of the differential equations for the case 
where the entrance Mach number is 2.8 and the tube wall is insulated were obtained with 
the new program using CS. The trial-and-error procedure set-up to find the correct 
Initial conditions la satisfactory. 

Programs are being prepared for the solution of the differential equations for 
the case where temperature dependence of the fluid viscosity and thermal conductivity is 
taken into consideration. 

T. Y. Toong 

Mechanical Engineering 
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-’04 N. EXCHANGE INTEGRALS BETWEEN REAL SLATER ORBITALS 

As is usually reported the testing of this program continues. However, with 
♦.he Insertion of an adequate set of constants for a sufficiently accurate guaaslan 
.quadrature, the program has proved able to give satisfactory results. Thus, at this 
juncture,the program has fulfilled the projected specifications both as to mathematical 
procedure and programming accuracy. The testing that shall be continued will be prima¬ 
rily devoted to more extensive Investigation ol' the complicated mathematical development, 
mat Is. of an order of magnitude larger than one amenable to hand computation. 

The program that has been developed wa3 envisioned as a step in the calculation 
ol molecular wave functions, and while meaningful alone, has Its primary use In this con¬ 
text, For this reason the projected succeeding step is a program designed to compute all 
of the integrals necessary to the computation of a molecule and organized to take careful 
account of the "computing weights' of each operation. This program should be forthcoming 
within the next year and will be described more fully In a later report. 

P. Merryman 

Laboratory of Molecular Structure 
and Spectra 
University of Chicago 






2If- C. ULTRASONIC DELAY LINES 

This problem is concerned with maximization of the reflection path aperture of 
a polygonal ultrasonic delay line. 

The aperture of such a delay line is aetermined by the difference between the 
least ordinate of one set of points in a plane and the greatest ordinate of another set of 
points. Expressed symbolically we have two sets of ordinates 
T T T B B B 

(H 0 , Hj . H f , ), and (H q . H t , ..., H p ). 


The aperture then la given by A - min hJ - max H*. The H i , and hence, also 
A, are functions of the size, shape, and position of the polygon in a coordinate system. 

A, however, is independent of the position. The size, shape, and position are determined 
by the distances, R., of the facets from the origin and the angles which they make with 
the x-axis. These J variables are given initial values, and then any new values are 
determined by changes,ARj and </j, from the old ones. 

T B 

The equations for H, and H. are linear in the variables dRj and are algebraic 
expressions of trigonometric functions of the af <. The latter are quite linear over the 
amall range for the J, with which we are concerned, so that the approximate linear equa¬ 
tions can be used to good advantage. We thus arrive at & set of equations 


n T n " 1 

H, 9-a,,4R.f 


J-o 

n-1 


1J 


n-1 

J-o 

n-1 


l.t 


b o B y 

H, b H, + £— 
1 J=o 


b, ,4H, +■ ZL b', , v, 
J 1 ® 


, 1=0.p 


l-o,... ,p 




where n = number of 8ide3 of the polygon. The matrices of coefficients 

{.b^ , £b' lt | are computed by a program on Whirlwind I as the initial step of 

the problem. 


Since we want to maximize A with respect to only the shape of the polygon, we 
must impose other conditions on the variables. Also, certain physical requirements are 
m*oe of the angles. Altogether, these lead to four other equations which the variables 
must satisfy ana, in audition, IT the polygon la symmetric, there are $n more conditions 
For each additional condition one variable is eliminated. This is done on Whirlwind an 
the second step of the problem, in addition to a transposition of the matrix and the 
addition of four variables needed in the linear program. 
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and Dr Edward N. Lorenz, Department of Meteorology• 


W. D. Sellers 
Meteorology 


1 i C. TRANSIENT TEMPERATURE OF A BOX-TYPE BEAM 

T _ < 5 nmmaru Reoort No 39 a program for computing the transient temperature and 

used h t^obtaln transient temperature and stress responses for the samehes^reauUs^are to 
exposed to various Intensities or thermal radiation on one surface These “”_ t0 

ce correlated with experimental work recently conducted by Mr. S. J. Engel Of me Aero 
elastic and Structures Research Laboratory. 


Luclen A. Schmlt 
Aeroelastlc and structures 
Research laboratory 


1*3 L. EIGENVALUE PROBLEM FOR PROPAGATION OF ELECTROMAGNETIC WAVES 

The previous report on this problem applies at the present time. A set of 
practical results which can be compared with the experimental results will be obtained In 
me near future. 

H. B. Dwight 
Lincoln Laboratory 


luL B.N. AN AUGMENTED PLANE WAVE METHOD AS APPLIED TO SODIUM 

Mr. M. Saffran of the Solid State and Molecular Theory Group continues to test 
the programs mentioned in the previous Summary Report. 

It h 03 been found that the matrix element generation routine which has been 
used Is overly slow and largely Inaccurate; therefore, a new routine Is being written. 

M. M. Saffren 

Solid State and Molecular Theory 
Oroup 


149 N. LAMINAR BOUNDARY LAYER OF A STEADY, COMPRESSIBLE FLOW IN THE ENTRANCE REGION OF A 
TUBE 

In oonneotlon with the research on heat transfer to a stream of air flowing at 
supersonic speeds In a round tube, a theoretical Investigation of the characteristics of 
the laminar boundary layer In the entrance region of the tube has been carried out. The 
boundary layer equations of continuity, momentum and energy are to be solved for specific 
entrance Mach numbers and thermal conditions at the tube wall. 

Gill's method is used In the numerical solution of these equations. 

Solutions of the first three sets of the differential equations for the case 
where the entrance Mach number la 2.8 and the tube wall le Insulated were obtained with 
the new program using CS. The trlal-and-error procedure set-up to find the correct 
Initial conditions la satisfactory. 

Programs are being prepared Tor the solution of the differential equations for 
the case where temperature dependence of the fluid viscosity and thermal conductivity la 
taken Into consideration. 

T. Y. Toong 
Mechanical Engineering 
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JO** N. EXCHANGE INTEGRALS BETVEEN REAL SLATER ORBITALS 

Aa la usually reported the teatlng of thla program contlnuea. However, with 
the inocrtlon of an adequate act of conatanta for a aufflclently accurate guaaalan 
; jadrature, the program haa proved able to give satisfactory reaulte. T^us, at thla 
:tincture, the program haa fulfilled the projected apeclflcationa both aa to mathematical 
procedure and programming accuracy. The testing that ahall be continued will be prima¬ 
rily devoted to more extensive Investigation of the complicated mathematical development, 
that la, of an order of magnitude larger than one amenable to hand computation. 

The program that haa been developed wa3 envisioned aa a atep In the calculation 
of molecular wave functions, and while meaningful alone, haa Its primary use In thla con¬ 
text For thla reason the projected succeeding step is a program designed to compute all 
of the Integrals necessary to the computation of a molecule and organized to take careful 
account of the computing weights of each operation. Thla program should be forthcoming 
within the next year and will be described more fully In a later report. 

P. Merryman 

Laboratory of Molecular Structure 
and Spectra 
University of Chicago 


216 C. ULTRASONIC DELAY LINES 

Thla problem la concerned with maximization of the reflection path aperture of 
a polygonal ultrasonic delay line. 

The aperture or such a delay line Is determined by the difference between the 
.east ordinate of one set of points In a plane and the greatest ordinate of another set of 
points. Expressed symbolically we have two sets of ordinates 

(H 0 T , H t T .H p T ). and (H^, H, B .H^). 


T R 

The aperture then is given by A • rain H t - max « 1 . The H 1 , and hence, also 
A, are functions of the size, shape, and position of tne polygon In a coordinate system. 

A, however, Is independent of the position. The size, shape, and position are determines 
by the distances, R., of the facets from the origin and the angles which they make with 
the x-axls. These J varlables are given Initial values, and then any new values are 
determined by changes. ARj and 6 j, from the old one3. 

T B 

The equations for H^ and H. are linear in the variables 4Ri and are algebraic 
expressions of trigonometric functlons^of the^,. The latter are quite linear over the 
small range for the J, with which we are concerned, ao that the approximate linear equa¬ 
tions can be used to good advantage. We thus arrive at a set of equations 


T ^ T n-2 

n-1 




Hi = \ 

j:o 

a H AR.t racL 

J J-o 

a 'lj 

, 1=0,. . 

. *p 

n-1 

n-1 




a O B y 

H. = H, C— 

1 J=° 

b. , AK. + ZL 

J J=o 

b 'u 

, 1-0,.. 

•. ,p 


where n = number of sides of the polygon. The matrices of coefficients , 

{a {.bj ] , [b' 1; l are computed by a program on Whirlwind I as the Initial step of 

the problem. 


Since we want to maximize A with respect to only the shape of the polygon, we 
must Impose other conditions on the variables. Also, certain physical requirements are 
made of the angles. Altogether, these lead to four other equations which the variables 
must satisfy and, in audition. If the polygon la symmetric, there are $n more conditions 
For each additional condition one variable is eliminated. This Is done on Whirlwind as 
the second step of the problem, In addition to a transposition of the matrix and the 
addition of four variables needed In the linear program. 
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H. D. Sellers 
Meteorology 


V,y C. TRANSIENT TEMPERATURE OP A BOX-TYPE BEAM 

Reoort No 14 a program Tor computing the transient temperature and 

MMkSpiimr- 

be correlated with experimental work recently conducted by Mr. S. J. Engel oi 
elastic and Structures Research Laboratory. 

Luclen A. Schmlt 
Aeroelastlc and Structures 
Research Laboratory 


1*3 L. EIGENVALUE PROBLEM POR PROPAGATION OP ELECTROMAGNETIC WAVES 

The previous report on this problem applies at the present time. A set of 
practical results which can be compared with the experimental results will be obtained in 
the near future. 

H. B. Dwight 
Lincoln Laboratory 


1*U R,N. AN AUGMENTED PLANE WAVE METHOD AS APPLIED TO SODIUM 

Mr. M. Saffran of the Solid State and Molecular Theory Oroup continues to test 
the programs mentioned in the previous Summary Report. 

It has been found that the matrix element generation routine which has been 
used is overly alow and largely inaccurate; therefore, a new routine is being written. 

M. M. Saffren 

Solid State and Molecular Theory 
Group 


199 N. LAMINAR BOUNDARY LAYER OP A STEADY, COMPRESSIBLE PLOW IN THE ENTRANCE REGION OP A 
TUBE 

In connection with the research on heat transfer to a stream of air flowing at 
supersonic speeds In a round tube, a theoretical investigation of the characteristics of 
the laminar boundary layer In the entrance region of the tube has been carried out. The 
boundary layer equations of continuity, momentum and energy are to be solved for specific 
entrance Mach numbers and thermal conditions at the tube wall. 

Gill's method is used in the numerical solution of these equations. 

Solutions of the first three sets of the differential equations for the case 
where the entrance Mach number is 2.8 and the tube wall is insulated were obtained with 
the new program using CS. The trlal-and-error procedure set-up to find the correct 
Initial conditions is satisfactory. 

Programs are being prepared for the solution of the differential equations for 
the case where temperature dependence of the fluid viscosity and thermal conductivity la 
taken Into consideration. 
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. 0*» N. EXCHANGE INTEGRALS BETWEEN REAL SLATER ORBITALS 

As is usually reported the testing of this program continues. However, with 
the insertion of an adequate set of constants for a sufficiently accurate guasslan 
uadrature, the program has proved able to give 8atlafactory results. Thus,at this 
juncture, the program has fulfilled the projected specifications both as to mathematical 
, rocedure and programming accuracy. The testing that Bhall be continued will be prima¬ 
rily devoted to more extensive investigation of the complicated mathematical development, 
hat is, of an order of magnitude larger than one amenable to hand computation. 

The program that has been developed was envisioned as a step in the calculation 
o! molecular wave functions, and while meaningful alone, has its primary use in this con¬ 
text. Por this reason the projected succeeding step is a program designed to compute all 
f the Integrals necessary to the computation of a molecule and organized to take careful 
account of the "computing weights' of each operation. This program should be forthcoming 
within the next year and will be described more fully In a later report. 

P. Merryman 

Laboratory of Molecular Structure 
and Spectra 
University of Chicago 


216 C. ULTRASONIC DELAY LINES 

This problem is concerned with maximization of the reflection path aperture of 
a polygonal ultrasonic delay line. 

The aperture of such a delay line Is determined by the difference between the 
least ordinate of one set of points In a plane and the greatest ordinate of another set of 
points. Expressed symbolically we have two sets of ordinates 
'p T T B B B 

(H 0 . Hi .Hp ), and (H q , Hj.H p ). 


The aperture then is given by A * rain Hj 1 - max Hj B . The Hj, and hence, alao 
A. are functions of the size, shape, and position of the polygon In a coordinate system. 

A, however. Is independent of the position. The size, shape, and position are determined 
y *.ne distances, R., of the facets from the origin and the angles which they make with 
the x-axls. These J varlablea are given Initial values, and then any new values are 
determined by changes. ARj andrfj, from the old ones. 

T B 

The equations for H. and H. are linear in the variables 6Rj and are algebraic 
expressions of trigonometric functions^of they,. The latter are quite linear over the 
small range for the J. with which we are concerned, so that the approximate linear equa¬ 
tions can be used to good advantage. We thus arrive at a set of equations 


a 'u <1 


1-0.p 


B o B y 

Hj Z Hj * <L_ 




b 'u J i 


where n - number of sides of the polygon. The matrices of coefficients Ja^J, 

[bj ( ] , [b'jjl are computed by a program on Whirlwind I as the Initial step of 
the problem. 

Since we want to maximize A with respect to only the shape of the polygon, we 
must Impose other conditions on the variables. Also, certain physical requirements are 
made of tne angles. Altogether, these lead to four other equations which the variables 
must satisfy and, In addition, IT the polygon Is symmetric, there are $n more conditions. 
For each additional condition one variable Is eliminated. Ttils is done on Whirlwind as 
the second step of the problem. In addition to a transposition of the matrix and the 
addition of four variables needed in the linear program. 
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linear Inequalities, rather than equalities. Tor which we wish to maxlmlie the quantity 
£ i + £ a - T 

c \ -^ij U J V J - ° H i 1 p 




U U J 


+ £d' 


1J V J 


& -°H, 


B 


1 = 0 , 


We irust make rurther substitutions, since the simplex method of linear programming Is to 
be used .""and Sis metnod requires that the variables and the constant part of the Inequal¬ 
ities must be non-negative. The variables y x and y 2 are Introduced Tor each uj. “e 
substitute a new variable Uj . Uj+yy and tot each Vj we substitute a new variable 
Vi • v< + y 2 . To make the constants non-negative we merely add .5 to each of them, and 
now consider our problem to be the maximization of ^2 - 1 e A. Since we know that 

It will also cone out In 


H, T and -H, B * -.5 Initially, this accomplishes our purpose. 

our solution that -y^ s min Uj and -yg = min Vj so that Uj and Vj are non-nega lve. 


An additional program Is required after the linear program to reconvert the 
results back to the original varlableB. The output of the results may then be provided 
for and an Iteration of the entire process may be set up In case the cfj are large enough 
so that they no longer can be treated as linear variables. 


The second step of this process actually requires two programs, one for the 
symmetric case and one for the non-aymmetrl: case. 

Since the first program requires a considerable quantity of data which can be 
determined Trotr. much less information, another program Is to be used to compute the data, 
and also to change the parameters in the other programs to accomodate different polygons. 


The first of the programs mentioned above,as well as the program for the 
symmetric case in the second step,have been written and checked, as has the linear 
programming routine. The latter, however. Is under revision to take advantage of the new 
si instructions for the auxiliary drum. The program used in the fourth step has been 
written. Some work has been done on the program for the non-symmetric case in the 
second step which will be similar enough to the symmetric case as to give no trouble. 

The program for setting up data la not started but will not be complicated. 


R. Biahop 

Arenberg Laboratory 


218 N. TRANSPORTATION OP INTEGRALS FOR DIATOMIC MOLECULES 

Several revisions were made In subroutines developed under this problem In 
connection with the work on Problem No. 288. 

R. K. Nesbet 

Solid State and Molecular Theory 
Group 
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COMPARISON OP SIMPLEX AND RELAXATION METHODS IN LINEAR PROGRAMMING 


A program has been written to ao.ve the classical transportation problem by 
the stepping stone method. 

The transportation problem may be stated as follows: A company operates plants 
producing » commodity the l 1 * of which can supply S, units or the commodity. The company 
sells Its production to n customers, the Jth of which desires D, units of the commodity. 
The cost or aanufactu-lng and transporting a unit of the commodity from plant 1 to custom 
iLh iLh e “ lre *-? rind the numDer of unit" that should be shipped from 

minimum*"' ***" 0u8tol,,er " ln ° rtl * r to have the total cost of the operation be a 
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Thus the problem may be stated mathematically as 
minimize C = ? C X 


subject to the constraints 


U> 


1,J 1J 1J 


^ X 1J - D J 

(2) 

f" x lj : S 1 

(3) 

*lj- 0 

(»> 


This Is a special case of the general linear programming problem. 


A sot of Xij which satisfies equations 2, 3, and 4 la called a feasible 
solution to the problem. if the set also minimizes C, It Is an optimum feasible solution 
containing. In general, exactly m + n - 1 non-zero Xjj. 

The process of solution consists of 

1) generating a feasible solution having exactly m»n - 1 non-zero X^. 

2) finding a zero Xj, which. If allowed to be positive, would yield a 
decrease ln C. In the process of Increasing this X t j, one of the 
non-zero Xjj must go to zero ln order that equations 2 and 3 remain 
satisfied. Thus the new feasible solution will again have exactly 
m+n - 1 non-zero Xjj. 

Step 2 Is repeated until there are no zero X, , which can be changed so as to 
reduce C. The set of Xj . is then an optimum feasible J solution. 

The Whirlwind program will handle problems for m t 128, in t- n p 400, and 
mn i r OOO. The program la presently In the final debugging stage. 


In the next quarter debugging will be completed and the program will be put ln 
convenient form for users. Production runs will be performed with a number of sets of data. 

This work la being supervised by Associate Professor W. K. Llnvlll of the 
Electrical Engineering Department, programmers are J. B. Dennis, Research Assistant, and 
W. J. Eccles, Teaching Assistant, ln the Department of Electrical Engineering. 

This project has been supported In part by a grant-ln-ald from the Union 
Carbide and Carbon Corporation. 


J. B. Dennis 
Electrical Engineering 


224 N. COMPUTATION 0? THE FIELDS OP VERTICAL VELOCITY AND HORIZONTAL DIVERGENCE 

The nature or the problem has been presented ln Project Whirlwind Summary 
Report No. 41. During this quarter the fields of vertical motion, horizontal divergence 
and vortlclty were computed for eight consecutive 12-hour periods. No additional compu¬ 
tations are planned until the Pressure Change Project of the Department or Meteorology 
has analyzed the results. 


J. M. Austin 
Meteorology 
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22b D. INVESTIGATION OP THE VOHTICITY FIELD IN THE GENERAL CIRCULATION OP THE 
ATMOSPHERE 

1. Description of Problem : The physical processes whlcn are Important In the maintenance 
of the general cTrculatIon of the atmosphere are being Investigated with the aid of a two 
parameter, non-linear, quasl-geostrophlc model or atmospheric flow. The model makes use 
of the vortlclty equation and the first law of thermodynamics and Incorporates effects of 
non-adlabatlc temperature changes and friction as well as the vertical advectlon of 
vortlclty and the transformation or horizontal vortlclty Into vertical vortlclty. A 
spherical coordinate system with pressure as the vertical coordinate has been used In the 
derivation, and the model Is being applied to the Northern Hemisphere rrom li*N. to 8CTN 
latitude. The two parameters used to describe the flow are the 700 mb and 300 mb contour 
heights. 

, „ . ^ th bbe a }? °f the »>dei described above It Is planned to examine the role of 

var.ous heating distributions In generating the gross features of the general clculatlon 
These distributions are described In the Initial report. (Quarterly Report Tor Project 

also £tven ?' Se °° nd Quarter - l‘A5). The equations of the model are 

also given In the initial report. 

- | r ? ce ? ure3 : An Averse matrix of a simple finite difference operator applied 
^ ^rld points has now been calculated. A solution involving the samp nimnia* 
has also been obtained by relaxation. The results are bllng fompareHn Je™s 

machine time used. It Is planned to repeat the work for 60 grid points and 

arour or non - hom °*™ — »«• 

aH3~The Sn® 0 "* 

Cambridge Research Center^y the^ol lowing 1 persone'^ der 3ponsorahl P of the *Ir Porce 

Richard L. Pfefrer 
Geophysics Research Directorate 

I>iant* S. Cooley 

Geophysics Research Directorate 
Paulo Castillo 

General Circulation Project, MIT 
Kirk Bryan, Jr. 

General Circulation Project, MIT 

Mr. Martin problem. 

D. S. Cooley 
Meteorology 


ATOMIC INTEGRALS 


whirlwind coding and applications 


A program ha. been constructed which evaluate. Integral, of the type 

IK U : A V* w . 

J74 3 N, Yu C.D . 

art Urary non-nega\ ive* Integer a ""a* and ^Tand'arbl trarv* ln Jff rals are evaluated Tor 
calculated wltr this program agiCe™^ ^d t eSi™iI*« P ° Bltlv# valuea of Results 

J. ‘•■?; 1 *i cant dsclBMils. Details or UiiJ prograS ill! bi° 2 fuIS ? 6t i een Bcven and •!*** 

Sol*d State and Molecular Theory Group. M.f.T., October*!^n^ln^Quarterly Progress Report, 

R. K. Nesbet 

Group State Moie cular Theory 

transient response op AIRCRAPT structure to aerodynamic heatino 

initiated o^Janunrv'^'inb’^bJVV'sctalt or'th^S*! £° aerod y n *»lc heating was 
Research Laboratory. A Program'toc^mput^L franlleA^ T' Aer ° ela ” tic «" d Structure 
butloi.fl lr: a structure simulating an airrr*rt w^ ran J len J tem Perature and stress dlatrl- 
The partial differential eauatloSa*2r hllf e? 1 8 haa been " rl tt«n and run successfully 
mattons similar to those repute?1 d / 8 ? J I ,d ** fl "“« <U™r™« 

realistic flight history was Incorporated n i®* 1 C!u ®: lerl y Progress Report. A more 

rlate and wet programs. incorporated In this program than In the thin plate or thin 

being planned? reo '‘ UD 0! ,r ' ls program are now being analyzed and several more runs are 

H. Parechanian 

Aeroelastlc and Structures Research 
Laboratory 

23'J C. 0UIDANCK AND CONTROL 

NO work was cone on this problem during the past quarter that was not classify. 

J. H. loaning, Jr. 

Instrumentation Laboratory 

241 B> "• TRAN3IENTS IN continuous distillation SYSTEMS 

ideal mouel or a contlnuoua*d 1stIllation coliitt" t ?hf tU S y i the transl *nt behavior of an 
has proceeded through the use ol thMedrSgiSZ; ThC 3tUdy ° f t "° 0lasaes ° r Problem 

:rs,“ “ sss.*}rt2 stsjs 

trdddet 118 ^T°°^’ ata ^ e ^“ a Pol : "l?d^d n «t?" art iff changed, and 

compositions are prevented from fluctuating ” !w"JlSlfie«t'S««.” lft * t lhe 

s. H. Davis, Jr. 

Chemical Engineering 
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c. DATA REDUCTION KOR X-1 EIRE CONTROL 

Thl« urotlem Is concerned with computing from gun fire control signals the 

“ islSrHSHw&w 

tt«» considered a termlnatlng report. 

The problem was coded by Dr. J. M. Slarx of the M I.T. Instrumentation Labora- 
tory. Results are c la a ail'lea Id format ion. 

.J. M. Stark 

Instrumentation Laboratory 


it N THEORY OF NEUTRON REACTIONS 

In the Summary Progress Report No. **2, a master program was described which for 
any given x will find the cross sections for any range of x‘ the U8e oi a auitabl<> 

parameter tape. ° 

Several logical errors have been uncovered ana eliminated in this master 
program. Tnere still remains a discrepancy between the cross-secttons computed by this 
: rogram and the known value Tor the particular case of u real square well. When this 
trour/le la cleared up It la planned to proceed as indicated In the previous progress 

report. 


E. Campbell 

Joint Computing Group 


N APPLICATION OP THE A t W METHOD TO BODY-CENTERED AND PACE-CENTERED IRON 

The energy band solutions for one-electron wave function In face-centered and 
i u.jy - e ntered Iron la proceeding via the Al’W (augmented plane wave) method. (Ref. 1.) 
Having chosen an approximate crystal potential, the first step is the formation of the bo- 
aiied Eve. Eo curves. (Ref. 1,2.) These curves are currently being produced for the 
fare-centered structure, using programs developed by Howarth and Saffren. 

A program is being written which will form the electronic charge density from 
the one-electron wave functions. 

J. H. Wood 

Solid State and Molecular Theory 
Oroup 

References 
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■' «1< C. WW1-1103 TRANSLATION PROGRAM 

A version of the two-pass Input translation program described In Summary Report 
No. 1 la operating satlsfactorily on WWI. Problems coded for the ERA 1103 computer In the 
described mnemonic code can now be fed Into WWI In the same automatic manner as CS flexo 
tapes with the rcultlng output being a machine-coded binary tape which can be read directly 
irit. an Ilf 1 -,. This procedure Is being used now by Problem No. 126 programmers, who are 
making use of both WWI an ; the Air Force Armament Center 1103 computer al Eglin Air Force 
Base. 

f°nv* modifications have been made to the translation program during the past 
quarter to make the binary tapes more easily acceptlble to the ERA photoelectric tape 
reader equipped llGVu Many more changes have been made to Improve the techniques used 


r ’rans Wit ion programs as well as to enlarge the vocabulary of the Input language, 
jumwtrv report on the program will be written In the next quarter. 

J. M. Frankovlch 

Digital Computer Laboratory 


YNAMIC ANALYSIS OP AN AIRCRAFT INTERCEPTOR 

During the past quarter, four solutions of thlB problem were run on Whirlwind 1. 
[' it. program consisted of a solution by the Gill method of eight nonlinear dlffer- 
, > •! • at Ions ueacrlblng the motion of an aircraft. The current runs solved thl3 basic 
; r . in, with certain changes in the original equations and In the Input constants and 
. t , u. The modifications were ascertained by a study of the results of runs made In 
i.i .?• , rev lows quarters. 

The .omplctlon of the present series of runs brings to an end th# Whirlwind 
• of the problem. Results of this phase and of the analogue study completed on the 
m :.t Flight Simulator of the Dynamic Analysis and Control Laboratory will be described 
PA L. Technical Retort No. 95, September, 1955. and W.A.D.C. Technical Report 
N -2*4', July, 1 . being written by engineers of the D.A.C. Laboratory and the 

•*l irt.ment of Aeronautical Engineering. 


K. Kavanagh 

Dynamic Analysis and Control 
Laboratory 


L. THE MEDIUM FREQUENCY IONOSPHERIC PROROGATION STUDY 

This problem was described In Summary Report No. 42. Data processing continued 
:ui'ir.g the pant quarter and la nearing completion. 

D. 0. Brennan 
Lincoln Laboratory 


f'U N. ELECTRONIC ENERGY OF THE OH MOLECULE 

Using the one- and two-electron lone- and two-center) Integrals, calculated on 
Whirlwind 1. the matrix elements lor the calculation of energy have been computed follow¬ 
ing the Lowain procedure for hand 1lng,non-oi-hogonal functions. These are now being UBed 
to calculate the ground state energy of the OH molecule at the observed lnternuclear 
distance. The 10 x 10 secular equation Is diagonalized on Whirlwind. 

A whole new set of basic integrals la now being calculated al an lnternuclear 
distance larger than the observed one. This will serve as a check on the previous calcu¬ 
lation and will determine the energy minimum of the E vs. r curve for the molecule. 

A. J Freeman 

Solid State and Molecular Theory 

Group 

261 C. FOURIER SYNTHESIS FOR CRYSTAL STRUCTURES 

The problem, which was first mentioned and described 1n Summary Report No. 42, 
was continued during this quarter. Particularly, the complete i -ulmensIona j latterson 
function for Diglyctnhydroch. 1 -.ride waa computed by Whirlwind from which, at the moment, 

Uv- minimum function' Is being obtained. Further, the crystal structure of Bectollte 
waa essentially solved, ar*u progress was made towards the solution of Wollaotonlte. 

Professor M. J. Buerger 

Geology and Geophysics 


4 2 


APPROVED FOR PUBLIC RELEASE. CASE 06-1104. 




WHIRLWIND CODINO AMD APPLICATIONS 

262 c. EVALUATION OP TWO-CENTER MOLECULAR INTEGRALS 

The evaluation of two-center Integra labelveen la iggpii^nav^been evaluated by 

bltals le partly completed. Two hundred an numerical integration program. Alao, one 
the uae of P. J. CorbatA's ^-generation * n J,"“"“"?;ited byP. Merrymin'a exchange- 
hundred and twenty exchange lntegrala have been evaluated oy r. mr y 

Integral program. 

This problem was deaorlbed in Summary Report No. *11. 


H. A. Aghajanlan 

Solid State and Molecular 

Theory Qroup 


26u C. OPTIMIZATION OP ALTERNATOR CONTROL SYSTEM 

nroareaa reDOrt for this problem, It was suggested that the prob- 
, r ontimizina the design of a system could be phrased as a minimization problem. Part 
of t^ work^completed In the past Sorter con.l.teS In coding a mlnl«lx.tlon procedurt B 

srswa 

rlvatlvea la Impractical. Therefore, a program was coded and wcoeMrully run «hlc(i p«r- 
forms the method of steepest descent and uses difference techniques to evaluate the de- 
rlvatIves. 

A scheme hae been worked out to apply the steepest descent technique to the 
problem of minimizing a function of n variables subject to constraints. fading of this 
scheme has not started yet, but will oe undertaken during the next quarter. Personnel for 
thlB quarter were R. M. Saundera, Profeaaor of Electrical Engineering (visiting), and 
J. B. Dennis and R. P. Newc, Research Assistants In Electrical Engineering. 

The work reported here Is under the sponsorship or the U. S. Air Pores, 

Contract AP 33(616)266*, Study of Aircraft A-C Oeneratlng and Voltage Regulating Systems. 

J. B. Dennis 

Eleotrlcal Engineering 


265 L.- ELECTRON D1PPUSI0N IN AN ELECTROMAGNETIC FIELD 

Work on tnls problem waa described In some detail in Summary Report No. 62. 

On a subsequent set of parameter values It developed that a smaller interval 
size waa required. Partial results were obtained for this set of parameters and the prob¬ 
lem waa terminated. 

D. N. Arden 

Digital Computer Laboratory 

266 A. CALCULATIONS FOR THE MIT REACTOR 

The program to calculate the reactor behavior for slow but sustained changes 
In reactivity waa successfully finished. The method outlined In the last quarterly re¬ 
port wsa used. 

A great number of data were obtained for aubcrltloal and supercritical systems. 
A strong variation In the average neutron lifetime with e, the reactivity, and the amount 
of poison present, was found. 

TTie metnoa described In the last report to obtain a zero determinant worked 
only with the heavy-water-cooled reactors; In cases with light water the round-off errors 
became too large. Alternate methods are being Investigated. All programming and 
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coding was done by M. Trooat. 

Tnla program la being carried out Tor the Nuclear Reactor Project. 

M. Trooat 

Chemical Engineering 

267 8. NUMERICALLY CONTROLLED MILLING MACHINE TURBINE BLADE 

The airfoil portion of the Jet engine turbine bucket uaed In this study la de¬ 
scribed by the coordinates of 272 points, the upper and lower boundaries of eight unequally 
spaced but parallel cross-sections each being defined by seventeen points. For any oroaa- 
sectlon the value of z la constant. On a cross-section boundary x la the Independent 
variable. When projected on the xz plane, theae form a trapezoidal grid. The bases or 
the grid are the traces on the xz plane of the root and tip cross-sections of the airfoil. 
The sldea are the xz projections of the leading and trailing edges. 

The total x Increment for each of the given cross-aectlone was divided by 
seventy-six and the si between each adjacent pair of croes-sectlone was divided by four, 
the result being 2233 Intersections. The intermediate z values were calculated by hand. 

A program, using linear Interpolation, was written to solve for the x values on oroaa- 
sectlons defined by the Intermediate z values. 

Another program waa written to solve the value of y at each of the 2233 Inter¬ 
sections for both the concave and convex surfaces of the blade. The method used was 
Newton's Interpolation formula In terma of divided differences. At the conclusion of this 
program 2233 x values for points on the concave surface of tnc bucket form a data block 
near the end of drum storage. The corresponding values of y form another block which 
starts after the last x value and runs to the end of drum storage. The twonty-nlne z 
values are held In magnetic core memory aa program constants. 

A third program waa written which evaluates the normal to the concave surface 
at each of the above pomtB with the exception of thoee on the edges or the surface. 

Each normal la the cross product of two vectors found by taking the first differences of 
x, y, and z between two points on either aide of point P on the surface. One vector la In 
the crosa-sectlon plane. The other la In the longitudinal section plane. The program 
finds coordinates of P', a point on the normal to the surface at F such that the distance 
from P to F equals the radius of the ball end mill uaed In the machinery operation on the 
Numerically Controlled Milling Machine (NCMM). Included aa constants In this program ara 
six reed rates In terms of Inches per minute. These were calculated by hand, It being 
decided from tool load conditions to change the feed rate at the end of the first, ninth, 
forty-seventh, and fIfty-aeventh passes of the tool. At the conclusion of thin program 
there are In drum storage 2025 data blocks of four parameters each. These are feed rate 
In Inches per minute and the x, y, and z coordinates of P'. These blocks are In proper 
sequence so that the tool center, when passing through consecutive points, will generate 
the convex surface of a die block that la In contact with the concave surface of the 
bucket during forging. 

With the series of four parameter data blooka In drum storage, uae waa made of 
a program written by A. Siegel of the Digital Computer laboratory Staff. The data blocks 
In storage are automatically translated Into NCMM Instructions by this program and the 
results are recorded on magnetic tape. After checking out the programs for Interpolation 
and evaluation of points along the normals, a run was made In which Siegel's program was 
the final phase. The ounchou paper tape obtained from the magnetic tape output was given 
a trial run on the NCMM using a wooden mock-up of the die block. This NCMM program ran 
satisfactorily for seventy-three of the seventy-five tool passes programmed. Near the end 
of the seventy-fourth pass the program broke down. The osuse was round to be In the com¬ 
puter program which evaluatee points on the normals. There are several straight-forward 
methods of modifying the programs that prepare the data for 31egel's program so as to 
allow the NCMM program to be obtained In Its entirety after a single computer run. Since 
time was not available to iw.ke a revision, It was expedient to recalculate the laat two 
tool passes by resetting certain constants and counters In the program at fault. The re¬ 
sults of this short run were spliced onto the NCMM tape at the beginning of the seventy- 
fourth pass. The corrected tape was rerun on the NCMM and the airfoil surface waa cut on 
a steel block. The results obtained were, substantially, those anticipated. 




This problem hae been the subject of a bachelor'a thesis entitled Numerical 
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Control for Machining Warped Sufj. The 
g UBWU T ST ts thA W^grimef. r ft !' ttlSTRFaa and 


problem has 
Engineering 


been terminated and the thesis 
Administration. 


0. Bromfleld 

Business and Engineering 
Administration 


2fQ B. CRITICAL MASS CALCULATIONS FOR CYLINDRICAL GEOMETRY 

The two-group, two-region method mentioned in the previous Summary Report has 
been used to calculate critical masses for uranlum-alunlnun-neavy-water 
by a heavy water reflector. On comparison with Snell s published data on this system, 
the calculatea critical mass averaged about 10* lower than the experimental one. Most or 
this difference can be accounted for by the shortcomings of two-group theory. The neutrons 
slow down much faster in diffusion theory than is actually the case; therefore, a two- 
group model does not take into account all the fast leakage. 


A parametric study of the effects of the fast diffusion coefficient and the age 
upon critical mass was made. As expected, critical mass was a strong function of the age, 
increasing with Increasing age. The correct critical mass could not be obtained without 
increasing the age to a fantastic amount, however. Surprisingly, the fast diffusion co¬ 
efficient had almost no effect. This Is probably due to the strong thermallzing effect 
of heavy water the fast flux drops almost to zero very quickly after leaving the core. 


To obtain better agreement with experiment, a three-group approximation 
(mentioned In the last report; will be tried. Also, a numerical solution of the trans¬ 
port equation by the S D method will be attempted, with a large number of regions and 
groups. As cylindrical geometry presents difficult problems for this method, a solution 
will be tested first in spherical geometry, and then extended to cylindrical cases. 

The two-group, two-region program is now being modified to obtain flux plots. 
Prom thesc, by numerical integration and perturbation theory, lifetimes and void coeffi¬ 
cients will be calculated. 


J. R. Powell 
Chemical Engineering 


271 B. EVALUATION OP A BEAM SPLITTING TECHNIQUE 


This problem is concerned with a statistical estimation of target position In a 
radar system when the operating conditions are such that the echoes received from the tar¬ 
get are corrupted by random noise and microwave interferen:e. A brief description of the 
method is given in the previous Summary Report and in the S. M. thesis report entitled 
Maximum Llkllhood Estimation of Target Position in a Radar System With Gaussian Noise 
Distribution’ submitted to the Department of Electrical Engineering in August, 1955, by 
Peter P. Engel. 


The quantitative results of interest were the means and variances of the fre¬ 
quency distributions of the estimated positions, and the variation of these quantities 
with the radar system parameters. The means of the frequency distributions were biased 
away from the true target position by approximately two percent, or leas, of the radar 
antenna's beam width, and the distributions had dispersions about the mean position of 
approximately seven percent of the antenna's beam width. The effect of doubling the 
antenna n angular velocity, thereby halving the sample size (or the number or received 
echo pulses;, was to increase the azimuth dispersion by approximately twenty-five per- 

° r h th * are presented In the above-mentioned thesis 

report, and the problem has been terminated. 


P. P. Engel 

Electrical Engineering 


27*1 L. GENERAL HAYDIST SOLUTION 

A program which derives the coefficients of an elgth degree polynomial and then 
proceeds to find a real root of the polynomial has been written. The program has proved 
to be quite contumacious and is still at the troubleshooting stage. 

M. Rotenberg 

Digital Computer Laboratory 


2~ \ N ANALYSIS OP AIR SHOWER DATA 

In Summary Report No. ^2 there la a description of the main features of the 
data analysis program for the air shower experiment. That report describes how the pro¬ 
gram la divided into a data processing part, a timing analysis, and a pulse height analy¬ 
sis. 

During the past three months the major errors in the program have been cor¬ 
rected and some testing has been aone with an artificial air shower and a real one. 

The data processing and timing analysis work satisfactorily. The main diffi¬ 
culties have been with the pulse height analysis. In the pulse height analysis we fit the 
fund ten f(a, A , X, Y) to the data by the method of least squares. The variables are nor¬ 
malized so that they are dimensionless and are all of the same order of magnitude. We 
denote by g* the numerical value of an event ♦'^corded by the 1 th detector. For each g. 
there is an equation of the form 

fit*. 4. X. Y) * g!♦<*!. 1 ■ 1. 2.N eq 

where 

*!<*. X, Y) * f(«, a, X, Y. X A , y A ) 

The numbers x^, are constants associated with the 1 th detector; ^ represents the un¬ 

known error in the measured quantity and also the error Introduced by the poorness with 
wnich the empirical function f represents the true functional dependence of air showers. 

It is necessary to Include in order to make the equations compatible for any given set 
of values of the variables CL, 0, X, Y. There are N e( , equations (N eQ ^ *0 and since there 
is a cf for each equation, the number of unknowns la w e q + ^ so 4 that there are in¬ 
finitely many solutions. According to the method of least squares, the best solution is 
the one for which the following expression Is a minimum: 

N e q N e q 

y (*• x, y)=£L *1 <*i 2 * *l( f i - 

i*«l x 1«1 1 1 


where w^ is a weight which expresses the reliability of the observed value g^. 


position 
function 
(best in 
process, 
is given 


We represent the four variables by a vector r » (<L, &, X, Y) so that f it a 
vector In a four-dimensional cartealan coordinate system and y Is a scalar 
of r. Let£L 0 , 0 O , Xo, Yo be initial guesses of the best values of d, A, X, Y 
the sense described above). The present program minimizes y by an iterative 
If r n »( a n . A n , X n , Y n ) Is the position vector after the iteration, 


grad (f n ) 

r " +1 * Pn '*|gra& y- (? n )f 


n “ 0, 1, 2,... 


where X is an arbitrary number which is chosen to make a rea sonable chan ge I n r per Itera¬ 
tion. It is required that y (r n +i) \fr (^ n ) and tnat gnfo gFatl y(r n ) .» 0. 

There are other conditions which must be satisfied but they are not Important and will not 
be described here. The minimization routine used is LSR lsi-201-271 (with slight modifi¬ 
cation). 


The main difficulty with this Iterative method is that it takes too much time 
to minimize y. At present it takes about three minutes for N e q * 20. It is hoped that 


^6 
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the time can reduced to one to a fora more suitable for calcu- 

of f and gPS) y by changing the <»Plrlc«l fJJJ? if, routine tray be necessary and 
latlon on a digital computer. ^ 8 ,iescrlbe* here mill be abandoned for some 

other J?oe\d^;:! bl irLr*p^niTo‘ii V vLriSlf. JSSS -thod described In rePerence 3. 
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280 B CORRELATION FUNCTION 

As indicated In the initial report, this problem concerns the application of 
i::t meorj of multiple time series prediction as developed by Norbert Wiener. From the 
j’rt•t teal standpoint, the nwat Important concept of the theory Is that an implicit fac¬ 
torization for single aeries permits 8 " factorization” for multiple aerlea. 

During this quarter, moat of the time has been spent In perfecting transform 
i : >*grams and the program for single aeries factorization. Program work remaining ln- 
lu.lcs multiple summation of multiple series, a recursion formula, and the prediction 

and error formulae. 

It Is noped tnat research on this problem will continue in the next quarter. 

P. Hanna 
Meteorology 


27U N. MULTIPLE SCATTERING OF NAVES PROM A SPATIAL AHRAY OP SPHERICAL SCATTERERS 

During the third quarter troubleshooting has continued on this problem. Pro- 
duct Ion runs are expected to begin soon. 

M. Karakashlan 
Joint Computing Group 


N. AUGMENTED PLANE WAVE METHOD AS APPLIED TO CHROMIUM CRYSTAL 

Radial wave functions for a chromium crystal potential have been obtained for 
M A vaiues and a grid of energies. Prom these.,starting valuea for integration have 
been obtained at a distance from the origin. This step shortens the time required for 
obtaining needed in the E vs. E 0 routine. 

07 

M. Saffron 

Solid State and Molecular 
Theory Group 


HORIZONTAL STABILIZER MODES, SHAPES AND FREQUENCIES 


A program developed under this problem number was used to determine the charac¬ 
teristic values (A) and eigenvectors Jx J for matrix equations of the type A • x = \ • B • x 


L. Schmidt 

Aeroelaatlc and Structures 
Research Laboratory 


278 N. ENERQY LEVELS OF DIATOMIC MOLECULES (LiH) 

In the preceding Summary Report we described an investigation concerning the 
electronic energy of the lithium hydride molecule which is in progress. Extensive use la 
being made of Whirlwind facilities In evaluating the overlap, one-electron, and two- 
electron integrals for this molecule for a number of lnternuclear distances, and alao for 
solving the resulting secular equations by diagonallzatlon of the Hamiltonian and overlap 
matrices which are formed. 

Results have been obtained for the complete configuration Interaction problem 
at the observed lnternuclear distance in which we have used the Hartree-Fock Is, 2s lithium 
wave functions, and a normalized Is Slater Atomic Orbital for hydrogen. The result gave 
about one-half the observed binding energy, and was not appreciably different from the 
two-electron problem In which the Is shell is kept filled. Therefore, we have extended 
the two-electron configuration interaction problem to several other lnternuclear distances, 
and it has been possible to construct a Morse curve quite accurately through the calcu¬ 
lated points. This phase of the work will be extended to several additional lnternuclear 
distances. 

T^e two-electron configuration Interaction problem has been enlarged to Include 
the possibility of one or both electrons occupying a lithium 2p orbital. Since the 2p 
energy levels lie near the 2s levels, It Is expected that the value of the calculated 
binding energy should be somewhat Improved by this treatment, although no results have been 
obtained at the moment, 


A. M. Karo 

Solid State and Molecular 
Theory Group 


288 N. ATOMIC WAVE FUNCTIONS 

Calculations on the states of C, Si, and Ce have been concluded. Similar 
calculations have been carried out, through the stage of evaluation of coefficients In 
the approximate linear expansion of self-consistent orbitals, for Mn and Vbf* , both in 
states. 

R. K. Nesbet 

Solid State and Molecular 

Theory Oroup 

291 B. DYNAMIC BUCKLING 

This problem is concerned with the unstable motion of structural columns sub¬ 
jected to large, t '.me-vary ing compressive loading. This motion may be in the elastic or 
plastic strain range of the column. 

The elastic-plastic response of the column after the loads have been removed 
has been approximated by the use of a 12 degree-of-freedom lumped maaa system. A second 
order of accuracy polynomial approximation Integration procedure is used to integrate the 
equations A program has been written and two successful runs have been made. Two addi¬ 
tional runs are planned to complete this phase of the work. 

The elastic response during the loading haa been studied by expressing the de¬ 
flection of the column by a combination of Its normal modes. The Kutta-0111 procedure has 
been used to Integrate the resulting equations. This motion has been calculated for four 
basically different variations of applied load versus time. One additional program may be 
run to secure additional data of this type. 

The work deacrloeo above haa provided useful Information about elastic-plastic 
vibrations and dynamic buckling response. It la anticipated that the next several days 
will complete this study. 

R. E. Jones 

Civil Engineering 
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296 C. SYSTEM ANALYSIS 


This program has been collecting data for a classified project at the , . 

l.ialrumen tat Ion Laboratory. Numerically the problem Involved the solution of seven simul¬ 
taneous ordinary differential equations using the Gill routine. 


This phase of the analysis has been successfully completed, but the system Is 
still In the development stage. Therefore, we are keeping the problem number open antici¬ 
pating that future computations will be necessary. 


w. B. Kohl 

Instrumentation laboratory 


297 B. DIFFUSION BOUNDARY LAYER 

A suggested solution to the thermal problem at the boundaries of high speed 
surfaces la the injection of a coolant through a porous wall. As part of a doctoral 
thesis, the appropriate system of equations describing a compressible laminar boundary 
layer flow for a binary mixture have been developed to determine the effects of foreign 
coolants. Mass and thermal diffusion are Included. 

By means of Blaslus and Dorodnitzyn type transformations, the ByBtem of equa¬ 
tions reduces to three ordinary differential equations of the form (thermal diffusion 
neglected): 

(Ar)‘* rr" = 0 



°i" ♦ [ f ♦ ( a 1 )3 & c i' • 0 
T" ♦ tjT' ♦ t 0 T . t ra 


which describe continuity of momentum, mass, and energy, f, c^, and T represent velocity, 
fractional concentration by weight, and temperature, while A relates the products of 
density and viscosity for the mixtures and external stream fluid. The remaining para¬ 
meters are functions of Prandtl, Schmidt, and Mach numbers, ao well as f, c i# the 
molecular weight ratio, and relative thermal capacities. Boundary conditions prescribe 
external stream conditions far from the surface, vanishing velocity at the surface, an In¬ 
jection rate, and, coupled with the latter, a concentration gradient such that no net flow 
of external fluid enters the surface. 

Several approximate treatments of the equations are possible. When the Schmidt 
number is identically unity, the concentration is linearly related to the velocity. In 
addition, when temperature effects upon the relative magnitudes of the thermodynamic 
transport coefficients are negligible, the momentum and mass continuity relations are In¬ 
dependent of the energy equation. Lastly, with an appropriate constant introduced Into 
the original transformation, one may assume A to be characterized by Its surface value and 
thus employ the profuse prior Information regarding the Blaslus equation. 

With the aid or Whirlwind, several of the mentioned simplifications are being 
evaluated using a Hellum-Alr mixture. Two point boundary conditions were replaced with 
Initial estimates at the surface modified successively so as to yield an array of in¬ 
fluence coefficients deduced Crom asymptotic results at Infinity, (infinity was essential¬ 
ly a value of 6.0 for the independent variable for asymptotic constancy to 6 significant 
figures.) The Glll-Kutta numerical integration procedure was employed. Programming In¬ 
cluded repeated Integration cycles using the more exact boundary conditions resulting from 
the influence coefficient system until end conditions were satisfied to 5 significant 
figures. Generally, convergence to the proper values occurred within four cycles. In 
some cases, the tabulated Blaslus function, f, and Its derivatives were uBed in conjunc¬ 
tion with a Taylor Series to yield the necessary intermediate values for the 0111-Kutta 
routine. 


The results indicate appreciable reductions in the recovery (thermal equilibrium) 
temperature for a flat plate at relatively low injection rates. Due to the presence of 
the T term in the energy equation, the recovery factor varies with Mach number. Ao a 
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result, for certain Injection rates and Mach numbers zero recovery factors are found. In 
all cases approximations to the exact equation yield conservative values. TTjIs is * con¬ 
sequence of the importance of the relative molecular weights for the components. Contem¬ 
plated work Is concerned with further examination of thermal diffusion effects as well as 
the limits for Injection rates. 

J. R. Baron 

Naval Supersonic Laboratory 


298 A. THE ELECTRONIC ENEROY OF THE HELIUM MOLECULAR ION 

The potential Interaction between a normal helium atom and a singly charged 
positive helium Ion is of considerable interest In connection with the understanding of 
the collision processes occurring In gaseous helium. This interaction has been considered 
by Pauling (1933) sod by Welnbaum (1^3b) who both used for the electronic wave functions 
of the * 2 y and states of He^ the symmetrized forms 


wnere X" is a doublet spin function given by 


X-(1, a. 

3) IS, )«j 

( 2 ) 


[IX -zr 


u(r|z) : 

4 

1 V 

(3) 


= u(r j |z* ) alrjjz 1 ) 

(<*) 


the suffixes A and B referring to the two helium nuclei and 1, 2, 3, to the three electrons 
of the system. 

Pauling took 2 - Z 1 and determined the electronic binding energy by the Rltz 
variational method, while Welnbaum carried the calculations one stage further by varying 
the two exponential parameters Independently. In the limit of large lnternuclear sepa¬ 
ration R it is^found that Z = 2 and Z 1 = 1.6875 which corresponds to vanishing interaction 
between the He Ion and the He atom. Iri the case of the ground state of the helium 
molecular ion Welnbaum found that at the equilibrium position the exponential parameters 
had the values Z = 2.02y and ■ 1.73**. which are not very different from those obtained 
in the limit of large R, and consquently, that the Improvement In the total electronic 
energy brought about by their variation was only slight. 

Welnbaum'a calculations correspond to using (*4) as the unperturbed helium wave 
function for large lnternuclear separation. In the present work (4) has been replaced by 
the expression 

^(r 1 ,r 2 ) - u (r 1 /a ) u(r 2 |,a ) + u(r 2 |a ) ) (5) 


where the parameters a. and fi have been fixed at the values 2.19 and 1.18 respectively, 
corresponding to infinite separation between the helium atom and Ion. This form for 
given a rather better binding energy for the helium atom than does expression (U) used by 
Welnbaum. As before, the total electronic energy for the molecular Ion can be determined 
by the application of the Rltz variational principle. The calculation Involves the 
evaluation of one-electron,coulomt . nybrld arid exchange molecular two-centre lrtegrala 
Tor a wide range of intemuclea- distances. These were performed by the Whirlwind elec¬ 
tronic computing machine using programs devised by Messrs. Corbstd, Dsigarno, and 
Merryman. The resultB of the calculation have not yet been fully analyzed, but will be 
available shortly. 


In conclusion. It Is a pleasure to express my thanks to Professor Morse, 
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299 C. HEAT TRANSFER IN TURBULENT FLOW 

The physical situation analysed wae that of a fluid flowing In well-developed 
turbulent flow lr. a smooth cylindrical pipe. The fluid initially at a uniform tempera¬ 
ture Tj enters a heated section having a constant wall temperature T w greater than Tj. 

It was desired to determine a relation between the heat transfer coefficient at any 
distance downstream and the physical and dynamic properties of the flowing fluid as charao 
terized oy the Prandtl and Reynolds numbers. 

In order to generalize the solution, tne following reduced variables are defined: 

fraction of possible temperature change 

T-T, 


0 - 


reduced radius 


and number of diameters downstream 

x = */D (3). 

The heat transfer coefficients are given in a generalized form in terms of Nusselt 
Numbers by the following expressions: 

. (R«) irr) ■ (4) 


Nu = ‘ 4 ( 1-8 )' 3T 1 


and the log mean Nusselt Number 


- ln(l-0 o > (5,) 


0 Q Is the nixing cup mean 8 at distance x given by 

f y u 8 dy 


where u Is the velocity at a reduced radlua y. The fractional temperature change at anv 
position Is given by the partial differential equation J 


$?] *> 
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with the boundary conditions 


B and W are functions of y only, as given by the following: 


in which u* Is a function of y* ( a function of y 

y* 1 (i-y) ^R|) Jt/z (14) 

The friction factor, f. Is a function of Reynolds number and may be 
approximated by 


f . 0.023(Re )‘* 2 ( 15 ) 

2 

The ratio or the eddy conductivity to the eddy viscosity, E/f , was assumed to 
The u*, y f relation used was as follows: 


be unity. 


y 21.995 
y + — 21.995 


u > = 13.56^5 tanh 




u r S 5.5 + 2.5 In y* - 15 

y* 


The constants were evaluated such that tne function and its Tlrst derivative 
are matched at a y* of 21.995. 

Solution of (7) was accomplished by means of a first order non-iterative 
rinite difference method. 

A variable spacing across the radius was employed for two reasons. One, the 
coefficients of the partial differential equation (7) approach infinity at the two boundaqr 
conditions, y = 0, and y = 1.0. Two, the so-called laminar layer at the wall Is very small, 
especially with high Reynolds numbers, and therefore small increments are necessary if the 
effect of the laminar layer la to be Included. If a stable solution is to be insured, very 
small x increments are also required, of the order ofay squared or less. 

In spite of the use of a variable y and x spacing, the number of repetitions 
required to obtain a solution is excessive. Ten tnousand repetitions of the calculations 
in the region of the smallest x Increment are required to calculate one diameter down¬ 
stream. This calculation would tax? about 15 hours on the computer. If a steady state Nu 
is to be determined, approximately 70 to 100 diameters distance down-stream must be calcu¬ 
lated . 

A 3 a test of the program, two very short runs lcomparer to the time required for 
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* complete Mint Ion) mere Me. One with the starting condition! epeolfled of 0 • 0; 
•m) another In which a temperature distribution was assumed similar to a condition some 
three quarters of the nay through the calculation. 


value. 

thesla 


The program ~ae successful. The results obtained were, or course, of little 
A complete description of the work In this problem will be presented In a Master 
In the rtepartment of Chemical Engineering. 


s 


A. Turano 

Chemical Engineering 


300 L. TROPOSPHERIC PROPAGATION 

A tape was cut for the correlation of radio signal strength versus layers In 
radio refractive index. A pilot program was completed, and further use will be made of 
the tape in the future for the analysis of other data. 

Values of a Fresnel integral have been computed, using a convergent power 
series with a complex argument. A number of setlefectory values have been obtained. 

H.B.Dwight 
Lincoln Laboratory 


3KU B. PARTIALLY CGWTDKJCUS WOCOW SEAMS 

For -ertair applicetlocs lr. the Department of Building Conatruction, Mr. M. 
Barker oes.rr.- .'inc the pro per: les of a wooden beam made up of two beams supported at 
the ends A and and M t-c .crier of overlays B. Hie two beams are connected at D and E. 
The overs., .engtr from A to C Is L. The other dimensions are fractions of L. 



‘ -» the o.stance :ros * tc the nearer end C, which varies rrom .25 to .50 In steps of .05. 

I o the end of a component beam at 0 or E, which Is to vary from 
•W He -‘-ope of .01. Y Is the maximu m deflection of the beam between A and B, 

wtint, .1 car be mtowr la the mutuum for the system. A uniform load W Is assumed. 

Pnr east taeae « cases, the forces acting at D and E. the bending moments, 
the shears anc T are to be found In terms of W, In which they are linear. 

—LI -eats to tee rational expressions In K and R of at most the eighth degree. 
The e.i.L-SM..,- Per T :s the maxim a of a fifth degree polynomial In the distance from A, 
whcee r:< 'T.c.erts oepecid on previously computed quantities. 


• the o«a-:.:.es that are used and those required have magnitude less then one, 
A-.C as a J .-r;..rv; is within the four decimal digit range of the Whlrxwlnu t „. 

ater. this ;rc:.ex au solve-; without recourse to the Interpreted CS computer. The 
;r Lmmrj rrsacr this prcc.em was presented to this Laboratory was to ssve time. A complete, 
sc :-rate ts.-.e o. About *50 numcers was returned to Mr. Barkan about s week after he 
^-vgmt^.- n.e lata and .ci»,a. It would have taxen him about a month using a desk 


303 B. PREDICTION OF CHROMATOGRAPHIC SEPARATIONS 

Very little Information la available for quantitative prediction of chromato¬ 
graphic separation by Ion exchange. Hawthorn performed experimental work for the rever¬ 
sible, non-equllbrlum case, where previous treatment la especially scanty. In addition to 
experimental work, Hawthorn developed a series of differential equations which describe 
the performance of an Ion exchange column effecting a c‘ romatographlc separation of two 
cations. Attempts to obtain solutions to these equatl. s were unsuccessful, but a method 
of numerical solution was presented. 

As an extension of Hawthorn's work, this thesis proposes the following objec¬ 
tives: 


1. Solution of tne equations derived by Hawthorn through the uae of the pro¬ 
posed numerical method and Whirlwind I. 

2. Comparison of these solutions, which have a theoretical baala, with the 
experimental data. 


3. Preparation of tables and graphs predicting elutrlation curves for chroma¬ 
tographic separations for various parameters. If objective number 2 shows good agreement. 

Re-examlnatlon of the basis Tor derivation of the differential equations. If 
objective number 2 shows poor agreement. 


In order to achieve the above objectives, a program of Instruction will be 
written for Whirlwind I. The first attempts will be kept as simple as possible, keeping 
close control on the adjustment of Increment slxe. etc. Due to the complexity of the 
equations, the numerical method used will be the Blmple Euler method where it Is assumed 
that 




or 


4u 

<Jx 


Is constant over the respective Increment. Convergence and stability of the solution will 
be checked by varying the Increment size. If the calculated values show good agreement 
with the nine sets of experimental data reported by Hawthorn, the Whirlwind program will 
be refined. If necessary, and production” runs will be performed for preparation of de¬ 
sign tables and graphs. 


The first two objectives listed above have been achieved. Solutions to the 
differential equations were round for two sets of parameters. Comparisons of these cal¬ 
culated results with experimental data give support to the validity of the assumptions 
used In the equation derivations. 

The following tapes have been left on file at the Digital Computer Laboratory: 

£ e \ Convergence test 

f« 303-232-20 ) 

fc 303232*86^ | Elutrlation curve prediction 


A complete description of the procedures used and the results obtained la 
given In the M. S. Thesis, "Prediction of Chromatographic Separations by Ion Exchange 
Using Whirlwind I" submitted to the Chemical Engineering Department In August, 1955, by 

J. L. Fischer. 


J. L. Fischer 
Chemical Engineering 


M. Weinstein 

Digital Computer Laboratory 


5 * 


55 

















APPROVED FOR 


WHIRLWIND CODING AND APPLICATIONS 


30b A. CALCULATIONS IN THE THEORY OP ATOMIC HYPERPINE STRUCTURE 

Measurements of atomic hyperflne structure (hfs) contain some of the bes * 
formation about the stationary electric and magnetic properties of atomic nuclei. In oroer 
to deduce accurate values for these nuclear properties (calleo nuclear moments!, one must 
have an accurate evaluation of certain Integrals over the electronic waveninctlons 01 tne 
atom, (Ref. 1). These Integrals have been evaluated previously (Ref. 2) by comparison with 
other atomic properties (the fine structure) using approximate but analytical wave- 
functions which neglect the screening of the nuclear electrostatic potential, we oecioeo 
to carry ou* a program of numerical computation of the wavefunctlons and Integrals, using 
a more realistic model for the atom. Since we require considerable accuracy In our 
results and expect to do the calculations for several atoms, tills undertaking appeared 
feasible only with the use of a high-speed computer. 


The equations to be solved are: 

where K, O. are constants; f(x),g(x) are the elgenrunctlons; £ Is the eigenvalue; and V(x) 
Is the potential seen by one electron In the atom. Por this potential we use (Ref. 3): 

v <‘> ■; [j^Txj 2 * *] 


where Z Is the atomic nucleus and ^ is t parameter which will be adjusted to give the 
"best" results. The variable x may range from 0to», but for physical and computational 
reasons, we have confined the study to .0001 i X S 12. This region Is divided Into b8k 
Intervals, with spaclngs varying from 0.000005 to 0.1. Numerical Integrations of the 
equations are done with the fourth-order Kutta-0111 subroutine. The ratio f/g at x = .0001 
Is determined from the boundary conditions that f, g be finite for x-»0. 

The main problem la to find the value of the eigenvalue ( which makes f and g 
go to zero as x goeo to •*> , Our method of finding e la patterned after that used by 
CombellCj (Ref. b). With some guessed value of < the equations are Integrated out to some 
large value of x, Xq (we use x Q = 12.). Por x a x 0 we set V = 0 and the solutions f, g 
for this outer region are then known In terms of spherical Hankel functions. The Inner 
and outer solutions are then Joined at xq: 

g (x 0 ), f(x 0 ) ~ Ci h (1) (lkx 0 ) -P Cg h (2) (lkx 0 ) 


where the order of the Hankel functions depends on the constant K and k rT 
Since the asymptotic value or these functions is: 

h^dkx)— e'^b 


h( z )(lkx) -» e kh 

the eigenvalue condition la Cg = 0. Thus a linear combination of f(x n ) and g(xo) multi¬ 
plied by Simple polynomials In (kx,,)- 1 gives the quantity Cge**0 whiofi we take as our 
error functions cf ( f ). Then repeating the Integration of the differential equations 
with new values or « we find the eigenvalue ( Q for which <f (£ 0 ) ■ o. with a fairly good 
Initial guess for < we can find £ 0 to about one part in 103 or better In 3 to 5 Iterations. 

It la probably worth noting that there are several ways In which one should not 
go about searching for the eigenvalue £ 0 . If one defines cf (£) from some such normalized" 
form as 

f(*o)/«(*o)> f(x 0 )/f(x 0 ), Cg/C lt etc., 

the denominator of each of the quotients has a zero very close to where passes through 
zero, and It Is very difficult to rind the zero of a function which looks like the 
drawing on tne next page. 
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Here the width of the Interesting region Is extremely small f*-' 10"^ In our problema). 

An alternative procedure used was to put In the desired (experimentally 
measured) value or £ and find the value of P (the potential parameter) that satisfies 
the eigenvalue criterion. When the eigenvalue le found, the several Integrals of 
Interest, 

/(r^g 2 ) dx, JTr 2 ** 2 ) x' 3 dx, ifgx' 2 dx, /rgx' 1 ‘dx, 

0 0 0 0 

are computed with the eigenfunctions using fourth-order integration. 

So far, theae calculations have been carried out for the p-doublet ground 
states of eight different atoms spanning the entire periodic table. The results are In 
some cases significantly different from the results of the earlier work. There are still 
a few more calculations to be carried out. A careful analysis of results In comparison 
with experimental data la expected to yield the values of the nuclear moments accurate to 
within a very small per cent. 
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305 A. COURSE 6.25. MACHINE-AIDED ANALYSIS 

Plfteen students In the summer section of 6.25, a course In maohlne-alded 
analysis in the department of Electrical Engineering, ueed TAC to solve a homework problem. 
The problem, Involving the generation of a sine function and the solution of a differ¬ 
ential equation, was designed to illustrate the use of an actual computer for solving 
numerical problems. As Is the usual case when students are first introduced to a digital 
computer, about 60$ of the students had programs work correctly the first time. 

W. J. Eccles 
Electrical Engineering 


306 D. SPECTRAL ANALYSIS OF ATMOSPHERIC DATA 

Recent theoretical, observational and experimental studies (conducted at M.I.T. 
and other Institutions) have demonstrated that the wave-like disturbances (eddies) In the 
atmosphere are of primary Importance In generating and maintaining the mean global rield of 
motion. To a large degree, the form of theae eddies characterizes the state or the general 
circulation. Little quantitative information has been obtained, however, as to the exact 
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scales or these atmospheric dlaturDancea and the coritrloutlon Which the airrerent scales 
or disturbances make to certain dynamically significant processes such as tne transport of 
heal arm momentum. 

Accordingly. it Is the purpose of this investigation to obtain long series of 
latitudinal space spectra for a number of Important atmospheric properties (“l™* 1 ® 
temperature, water vapor, vertical velocity, tranaport or neat and momentum) 88 a function 
of latitude, preanure-helght, and time. These spectra will serve as the basis for: 

1. Providing a picture of the time and space variations of the spectral 
distributions of significant atmospheric parameters, as well as a measure ol tne relative 
Importance of the various scales of eddies In that atmosphere. 

2. Testing the results of certain theoretical models of the atmospheric flow, 
regarding the stability properties of disturbances. 


1. Relating changes in the kinetic energy spectral distribution to other 
Important parameters. 

4 . studying the dynamic Interactions among disturbances of different scales 
and the mean flow. 


In more concrete terms, It is proposed to obtain the spectra for the 100U, 700, 
^OU, JUG, and ^00 mb pressure surfaces at latitudes 30°. and 60 N. 


Extensive computations of the type outlined above have become- icu m .u only 
w^th the advent of hLgh-apeed machines such us WWI. 


'IYie program is currently in the process of being checked. 


D. Saltzmari 

F. Was 

Meteorology 


307 C. SUPERSONIC NOZZLE DESIGN 

The Friedrichs method employed at the Naval Supersonic Laboratory Tor deducing 
supersonic nozzle contours has been programmed for Wt.iriwlnd*. On the ' asis of an assumed 
velocity distribution (nozzle generating function) along a reference axi3, there follow 
expressions for the physical coordinates and 1 i w jropertlea along streamlines for a 
perfect gas expanding from subsonic to supersonic sj•co It Is required that the exhaust 
fiuw be uniform and parallel for modei test purposes. Eacn dealgr. (exit) Mach number 
therefore necessitates a numerical integration process to establish the flow properties 
along a (design) characterl3tlc in the supersonic field. Downstream of the latter* charac¬ 
teristic line, streamlines are determined on a maos continuity basis such that a Prandtl- 
Weyer flow region precedes the final uniform flow field. 


The spectra will be based on data obtained at points (every degrees of 
longitude) around a latitude circle and will be computed for wave numbers 1 through 1' by 
means of the standard methods of harmonic analysis. In general, the complex spectral 
function for any property, q, Is given by 

P(n) = p,(n) + 1 p 2 (n) 


The Gill-Kutta procedure was employed for the numerical integration portion of 
the program. Application of the Nrwton-Rapnson method was made to find certain Mach num¬ 
ber distributions at equal coordinate Intervals. 

Tne program has been initially successful and Is currently oeing checked in de 
tall against earlier computations. It Is Intended tc make use of the program for future 
nozzle requirements at the Naval Supersonic Laboratory. 


where 


Pj(n) 



q( \ ) 00 s n A 


J. R. Baron 

Naval Supersonic Laboratory 


P 2 (n) 



1 


Here X represents longitude and n is the wave number. 


(A=o°, i°, io°,...3D5°) 


By suitably combining the Fourier coefficients for the heights of the pressure 
surfaces, temperature, the zonal wind, water vapor and vertical motion it is possible to 
obtain ail of the desired transport spectra and the kinetic energy spectra. For example, 
the spectral function, T(n), for the meridional transport of sensible heat for a given 
latitude pressure p, and time t, is given by 

T t,P,t< n > = n /M n ) a l< n ) - Mn) a 2 (n)] 'f.P.t 

where a, and a 2 are the real and imaginary parts, respectively, of the pressure-height 
spectrum, and b\ and b 2 are the components of the temperature spectra. 

The spectral function for the geostrophlc kinetic energy of the north-south 
flow is given simply by 


K f,p,t (n > Z n ‘ f a i< n r v a 2 (n) 2 J . 


It 

lylty, as well 
be processed. 


18 anticipated that complete hemispheric data for the month of January, 
as a much longer series of height data for the 500 mb level alone, will 


306 A. FREQUENCY ANALYSIS OP APERIODIC FUNCTIONS 

This problem is a continuation of problem lybC. We now plan to generate cosine 
functions of different frequencies ann cross correlate each function with the Intestinal 
motility data. We then will calculate the root mean square values of these cross correla¬ 
tion functions. 

J. Roseman 

Digital Computer Laboratory 


30b> B,N. PURE AND IMPURE KCI CRYSTAL 

A theoretical study is to be made of the electronic structure of an alkali 
halide crystal. The free-lon, Hartree-Fock wavefunctlons will serve as a starting point 
of the study. Since an alkali halide lr> Its ground state la a closed shell system, the 
ground state wavefunctlon may be approximated by a single Slater determinant of the nor¬ 
mally occupied, free-lon functions. Due to the small overlap of these functions In the 
crystal, the approximation Is Talrly good. If the determinant la made up of linear com¬ 
binations of these free-lon functions, the crystal energy is of course unchanged. For the 
perfect crystal It Is proposed to Determine one-electron wavefunctIona as those particular 
orthonormal linear combinations of the free-lon functions which satisfy the Bloch condi¬ 
tion and which satisfy the Hartree-Pock equations. From this the Hartree-Fock energy 
band structure of the crystal will be obtained and can be compared to recent X-ray studies. 

It 18 further proposed tu study the V-cent*.*r no a particular example of the Im¬ 
pure or almost perfect crystal. The V-cehter in an alkali halide crystal is a positive 
Ion vacancy with an associated, bound hole (an electron vacancy in an otherwise closed 
shell structure). The many-electror. wavefunctlon can again be approximated by a single 
determinant which is made up of linear combinations of free-lon functions. Using an 
argument due to Koopmans, however, this wavefunctlon *.a completely determined when the 
coefficients of only the missing electron are known. Thus the crystal energy for various 
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Staten of excitation of tne hole can be determined by solution of a finite secular equa¬ 
tion. Because of trie localized nature of the Impurity, the secular equation need only 
involve free-ion functions on sites near the Impurity center. This calculation does not 
depend on any preliminary, perfect crystal calculation. It Is thus proposed to determine 
transition energies for the bound hole and compare the results to recent experimental 
observations. 

In both the pure ana impure crystal problems described above, we must solve a 
set of simultaneous eigenvalue equations of the following form: 

nf C H nm“^Aim) a m * °* 

where the a m ‘s are the coefficients of free-ion functions or linear, symmetry combinations 
of free-ion functions, and H nm and A nn , are respectively energy and overlap Intervals. The 
Integrals are closely related In tne two problems of this paper; hence the problems are 
being performed together. 

Many of the Integrals Involved are oetween functions on different centers In 
the lattice. These Integrals are to be evaluated by the expansion method of Barnett and 
Coulson as adapted for Whirlwind use by P. J. Corbatd of this group. In this method, an 
Integral between wavefunctlons on different centers Is calculated by choosing one of the 
centers and expanding all of the displaced functions In a series of radial functions 
(cal led<l(r) by Corbatd) times spherical harmonics about that center. For two-electron 
Integrals l/rjg I s also expanded In spherical harmonics about the common center. The 
original Integral then becomes a finite or Infinite sum of integrals, each of which can be 
calculated easily. 

The order of calculation Is as follows. The Hartree-Pock radial functions for 
Cl* and K‘ are fitted with sums of exponentials times powers of r. The resultant analytic 
functions are expanded about another center using Corbatd's Alpha-Generation Program 
(named for the coefficients a(r)). Certain linear combinations of the alphas are per¬ 
formed by my own Mode Addition Program. Finally the basic integrals are calculated by 
Whirlwind using Corbatd’a Integrator Program. The basic Integrals are then combined by 
hand to form the H™, andA__ mentioned above. Finally the secular equations for the a-'s 
and the *'s are solved on Whirlwind using routines developed by Meckler and Corbatd. All 
of the programs mentioned have now been tested and have been found satisfactory. It is 
therefore expected that by the end of the next quarter most of the production runs will 
have been completed. 

L. P. Howland 

Solid State and Molecular 

Theory Qroup 


310 C. TRAJECTORY CALCULATIONS POR A ROCKET DURING POWERED FLIGHT 

A rocket obeys Kepler's Laws when It Is In Tree flight, since no exterior 
forces are applied however, when the rocket's engines are providing thrust Kepler's Laws 
no longer apply. Consequently, In order to find the optimum powered-flight trajectory, 
l.e.. the trajectory which uses the least amount of fuel. It Is necessary to do a trial 
and error trajectory study. A further complication is the effect of the atmosphere on 
the rocket, namely the forces or drag and lift. The equations of motion of a rocket 
during powered flight, taking into consideration all the forces that have been mentioned. 

hpp r.hf*n 


ip n (it) z - % *•'»/ + 


where the equations are written with reference to axes aligned with 
(Figure 1 defines the nomenclature uaed In the above equations.) 


the local vertical. 
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For the purposes of this Investigation the equations of notion were written 
in the following non-dimensional form 
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it typical trajectory la shown In Figure 2. In order to optimize the trajectory, 
various parameters must be varied, such as the length of time spent In vertical flight, 
the Impulse angle-of-attack at tip-over, the lerigtn of time before staging, the length of 
time before burnout, and the pitch-program used. Por any given set of parameters the 



velocity, flight path angle, altitude, and distance from launching point at burnout must 
be those which will cause the following equation to become zero. 
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311 N. SOLITARY WAVE GENERATING CAM 

This 1s the first quarterly progress report on Problem 311 N. The problem 
consists in the solution to a series of algebraic equations which describe the loci 01 
the tool centers for work to be executed by the numerically-controlled milling machine 
of the Servomechanisms Laboratory. The work for the milling machine la the cutting of 
stainless-steel cams to be used with the wave machine of the Hydrodynamics Laboratory, as 
part of the work carried out on the Solitary Wave Research Project, DIC 3-6ot>^, conducted 
for the Office of Naval Research under Contract Nc. Nt>orl-07837. 

The form of the cams la two half-length sine waves superimposed on a neutral 
circle of the cam. The ratio of the two wave lengths varies from 1:16 to 1:1 for a series 
or eight cams. The equations for the cam surfaces were modified to account for the center 
position of the milling machine tool during cutting. 

The input to the milling machine is binary punched tape; therefore, a Whirlwind 
subroutine developed by John Runyon was used to convert the decimal results to the alge¬ 
braic equations to binary punched tapes. 

The progranfners for this problem are Donald C. Taylor and John G. Houaley of 
the Hydrodynamics Laboratory. Technical assistance In programming was given by Mr. A. 
Siegel of the Digital Computer Laboratory staff. 

J. 0. HouBley 
Hydrodynamics Laboratory 


312 L. ERROR ANALYSIS 

The errors In the knowledge of a set of six parameters descriptive of a system 
are to be determined from the errors in certain measurements on the system. These measure¬ 
ments are assumed to be corrupted by additive gausslan noise, A small error approximation 
is made so that errors In the parameters are linear in the measurement errors. The para¬ 
meter errors are, therefore, also gausslan. Whirlwind 1 will compute the moment matrices 
of the errors In the parameters for various input conditions. 

The problem can be stated mathematically as the evaluation of a *4 x 6 and a 
6x6 matrix. The elements of these matrices are functions of six parameters and a 
variable of Integration, whose values will be determined from the limits of integration 
to be described below. The functions require the evaluation of square roots, trigonometric 
and inverse trigonometric functions. The parameters will be given as input conditions and 
the results determined for various sets of these parameters. The two matrices are to be 
multiplied together to yield a third matrix, the elements of which will be multiplied 
together two at a time and each pair integrated. The limits of Integration will be sup¬ 
plied and the results again determined for various sets of values. The integrals will 
form the elements of four 6x6 matrices which are to be added together with weighting 
factors which will also assume various values. The resulting 6x6 matrix la to be in¬ 
verted and then premultiplied by a 2 x 6 matrix and postmultlplled by the transpose of 
this matrix. The elements of the 2x6 matrix will be functions of the six parameters 
mentioned above and will require evaluation of the same basic functions. The final 2x2 
matrix is to be diagonalized. The Inverted 6x6 matrix, the eigenvalues and one eigen¬ 
vector of the 2x2 matrix are desired as final results. The various combinations of 
input values needed will require approxlmatel'' 100 repetitions for the basic procedures 
described. 

The program has already been tested and will be used in the near future. 


313 D. COURSE 6.601, NUMERICAL CONTROL OP MACHINE TOOLS 

As part of Course 6.601, the students are required to produce, using manual 
computation, a program for the Numerically-Controlled Milling Machine (NCMM). Because 
of time limitations, they are not required to carry their computations to completion but 
are permitted to stop at an Intermediate point at which the remaining calculations are 
obvious. 

A Whirlwind I routine was prepared which completed the students' computations 
and punched the resulting NCMM control tape. This routine was used to process all the 
data provided by the students, and the tapes were then used on the NCMM. The routine 
operated correctly, and produced satisfactory NCMM tapes. 

A tour of the WWI computer, during which was demonstrated the NCMM translation 
program developed under problem 2^0, was arranged as part of the course. 

A. Siegel 

Digital Computer Laboratory 


31^ C. FACT0RIN0 HIGH ORDER POLYNOMIALS 

Tftls problem concerns the analysis of the effects of aircraft flexure coupling 
with aircraft automatic control systems. The analysis results in two sets of six simul¬ 
taneous differential equations of motion. Whirlwind was used to determine the roots of 
the resulting performance functions. Factorization of several polynomials of about the 
10th degree was done using Hitchcock's method, for which programs were written by M. 
Jacobs of the Digital Computer Laboratory. 

V. W. Howard 

Aeroelastlc and Structures 

Research Laboratory 


316 L. RADAR CORRELATION 

A method la sought to correlate on a real time electronic computer the data 
that may be found directly by radar apparatus. 

Some time was spent to determine the analytic solutions which would have the 
least Inherent error and would allow high-speed computation. 

! To minimize error in a function of more than one variable, a linear programming 
routine is being developed. This routine has already been checked In part and is now 
being modified to take advantage of new instructions pertaining to access to the magnetic 
drum storage. By evaluating the exact expression at several points and then letting the 

I coefficients in a simplified expression be variable, we should be able to find coefficients 

that minimize the maximum error. 

M. Weinstein 

Digital Computer Laboratory 


I. I. Shapiro 
Lincoln Laboratory 
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NOTE: 

Reports on the following problems, bone by members of the Machine Methods of 
Computation group, may be found in Section 2.2 of Part 1. 


122 N. 

172 B.N. 
21 N. 
235 B.N. 
2U2 N. 


COULOMB WAVE FUNCTIONS 

ENERGY BANDS IN ORAPHITE 

ATOMIC WAVE FUNCTIONS AND ENERGIES 

EIGENVALUES IN A SPHEROIDAL WELL 

CALCULATION OF NUMBERS OF STRUCTURES OP RELATIONS 
ON FINITE SETS 


A. Terakin 
Z. Pried 
A. Tubla 

F. Corbatd 

A. Tubls 

J. Uretaky 

D. Mcllroy 
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APPENDIX 

1. SYSTEMS ENOINEERINO 


Magnetic Dr ums 

The Interlace on the auxiliary section or the buffer drum has been changed to 
provide more rapid access to drum storage. The computer may now read or record blocks of 
Information at the rate of one word every 32 microseconds. The former rate was one word 
every 6*4 microseconds. 

The control systems for both drums have been modified to permit the selection 
of a new mode of operation. The programmer, by adding 1000 (octal) to the present drum 
orders, may select the drum as a consecutively addressed storage medium with no discon¬ 
tinuities at the end of each group. The Group Selection Register (GSR) has been made a 
counter which la added to by the end carry from the Storage Address Register (SAH). The 
SAR end carry to GSR Is gated by digit 6 of the ln-out switch. 

Magnetic Tape 

In several Instances, valuable information recorded on magnetic tape has been 
destroyed by accidental recordings over this Information. A system has been Installed to 
allow the three printout units to be locked In the read mode by meana oi toggle switches 
located on their control panels. The computer may sense an Intervention bit to determine 
whether these switches have been thrown or whether the units are r.ady to record. 

Power Supplies 

A spare filament alternator has been obtained and the Installation work neces¬ 
sary for connecting It into the system Is in progress. The control system will be such 
that the alternator may be used to replace either of the two alternators now In service 
by operating one or two toggle switches. A system for substituting a motor-generator 
set for one of the d-c supplies la being developed. 

Polaroid Land Camera 


A Crown-Graphic camera with a Polaroid Land-camera back has been mounted In the 
control room for photographing the Indicator lights. The photographs may be used In place 
of manual recordings of the lights. This system provides a sure, fast, accurate means of 
gathering trouble-location information. 


Computer Reliability 17 June 195b to 22 September 1W55 


Total computer operating time 

2017 

hours 

Total lost time 

5*4.1* 

hours 

Percentage operating time usable 

97. *• 

percent 

Average uninterrupted operating time between failures 

1*4.2 

hours 

Failure incidents per 2**-hour day 

1.65 


Average lost time per incident 

23.6 

minutes 

Average preventive maintenance time per day 

1.5 

hours 










APPROVED FOR PUB 


2. VACUUM TUBES 

Vacuum Tube Life 

During the first three quarters of 195&* the Whirlwind computer operated Tor 
1/tiO, 1932, and 1800 hours, respectively. Vacuum tube life has been calculated for the 
six most numerous types In WWI, comprising approximately 600 of the total tube complement. 
A summary of this Information is shown below: 




FAILURE 

HATE. PER CENT 

1000 HOURS 


Tube Type 

Quantity 
now in 
Service 

1904 

Average 

Plrst 
Quarter 
£22_ 

Second 

Quarter 

_ 

Third 

Quarter 

_ 

7AD7/SR1407/6145 

4273 

1.5 

0.97 

1.1 

0.69 

7AK7 

3400 

0.44 

0.29 

0.18 

0.08 

5965 

887 

0.32 

0.38 

0.52 

0.60 

6080/ 6080W A.' 6 AS7 0 

712 

1.1 

0.48 

2.0 

3.5 

6BL70T 

496 

0.41 

0.46 

0.52 

0.96 


The number of 2D21 failures nas been materially reduced as a result of circuit 
changes correcting excessively negative grid voltage during conduction (see Summary 
Report No. 40, Pourth Quarter, 1904). 

The number of 5963 failures has been steadily decreasing. This Is because most 
of the older tubes of the original complement have already failed due to cathode inter¬ 
face Impedance and have been replaced; the present failures either went into service aTter 
the original equipment installation or are in circuits where operating margins are wider 
ana hence more tolerable of tube deterioration. 

The type 5965 continues to be quite stable, from a maintenance-of-characterlstics 
point of view practically all the failures In this ty^e are due to high negative grid 
current, primarily leakage, with some grid emission beginning in the older tubes. 

The relatively large fluctuations in the failure rate of 7AK7'8 are caused by 
statistical variations In the small number of failures per quarter. The failure rate for 
this type continues to be comfortably low. 

The type 6BL70T has now reacned approximately 12,000 hours of operation. Although 
the failure rate has been relatively low, it has been Increasing as a result of a small but 
lncreaalng number of low-plate-current failures partially due to cathode Interface Impe¬ 
dance. The sharp Increase In the failure rate during the third quarter of 1955 Is a result 
of the routine retesting of a panel temporarily removed from WWI service and the subse¬ 
quent failure, at test, of seven tubes for low plate current. The circuits in which 
these lubes operate (cathode follower being driven by a cathode follower) are such that a 
considerable decrease in zero-bias plate current can be tolerated before system failure 
begins. The extremely low zero-bias plate current cannot be completely ascribed to the 
magnitude of cathode interface Impedance involved (100-200 ohms). Qualitative tests In¬ 
dicate a major decrease In cathode emission Itself. 

Some comment Is In order regarding the high failure rate of the types 
6060/6080WA. A considerable quantity of these tubes Is used in series - tube - regulator 
panels where, prior to a recent modification, circuit conditions were such that the tubes 
were being operated considerably over plate dissipation rating. In the magnetic core 
memory system where operating conditions are much less severe, the failure rate of this 
type continues to be reasonably low. 
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3. PUBLICATIONS 


Project Whirlwind technical report* and memoranda arc routinely dlatrltuted 
to only a restricted group known to have a particular Interest in the Project and to 
ASTIA (Armed Services Technical Information Agency) Document Service Center, Knott Glid¬ 
ing, Dayton, Ohio. Requests for copies or Individual reports should be made to ASTIA. 

The following Is a 11s. of memoranda putllshed Ly the Scientific and Engineer¬ 
ing Computations Group during the past quarter. 


No. 





DCL-Sl 

Abstracts of memos on Automatic Coding 
Techniques Used on WWI 

6-22-55 

J M. 

Prankovlch 

DCL-82 

A Proposal for a New Method of Select¬ 
ing the WWI Drums for Block Transfers 
of Information 

6-21-55 

J. M. 

Prankovlch 

DCL -83 

Interpreted Instruction Code of the MIT 

CS II Computer 

6-24-55 

J. M. 

Prankovlch 

DCL-84 

The MIT Comprehensive System of Service 
Routines 

6-27-55 

P C. 

Helwlg 

DCL-88 

Normal S & EC Operating Procedures 

7-15-55 

P. C. 

Helwlg 

DCL-89 

Recent Changes In the CS Conversion 

Program 

7-15-55 

J. M. 

Prankovlch 

DCL-90 

Whirlwind Display Program: Vibrations 
in a Length of String 

6-17-55 

D. N. 
Q. W. 

Arden and 
Patterson 

DCL-91 

Modification of the Numerically- 
Controlled Milling Machine Translation 
Routine for Use with the G&L Machine 

7-19-55 

A. Siegel 

DCL-92 

Proposed Additional Function for the 

Read-In Button 

7-22-55 

F. C. 

Helwlg 

DCL- 94 

Polynomial Factorization 

8-1-55 

M. Jacobs 

DCL-94-1 

Revision In Polynomial Factorization 

8-17-55 

M. Jacobs 

DCL-97 

Whirlwind I Processing Routine for 

Course 6.601 

8-18-55 

A Siegel 

DCL-101 

Iterative Solution of Linear Systems 

Having Sparse Matrices 

9-12-55 

M. D. 

Mcllroy 
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U. VISITORS 

Tours of the WWI installation Include a showing of the film Making Electrons 
Count , a computer demonstration, and an Informal discussion of the major computer com¬ 
ponents. During the past quarter, 11 groups totalling 253 people visited the computer 
Installation. Included In these groups were: 


June ?8 


July 11 
August 3 
August 29 
August 31 


September 9 
September 13 
September 22 


Training Course for Whirlwind I 

MIT Porelgn Students Summer Project 

Course in "Machine-Aided Analysis" 

Westlnghouse Science Teachers Project 

Training Courae for Whirlwind 1 

Summer Session Course In "Numerical Control 
of Machine Tools." 

MIT Textile Department 

MIT Preshmen 

Worcester Polytechnic Institute Seniors 


The procedure of holding Open House at the Digital Computer Laboratory on the 
first Tuesday of each month haa continued during this period. Three groups totalling 35 
persons visited th- Laboratory at the Open House demonstrations. These persons represent¬ 
ed members and Trlends of the MIT community, Boston University, Harvard University, Tufts 
College, Royal Canadian Air Porcc, Michigan Department of Revenue, Arenberg Laboratory, 
and the University of Illinois. 

During the past quarter there were also ^3 Individuals who made brief tours of 
the computer Installation at different times. Represented y these Individuals were 
Reliance Electronics and Engineering Co., Caswell Aeronautical Laboratory, U. S. Army Ord¬ 
nance, Pratt and Whitney Aircraft Co.,Convair, and the Air Force Armament Center. 
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